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Abstract The a-particle condensate in nuclei is a novel state described by a prod- 
uct state of a's, all with their c.o.m. in the lowest 05 orbit. We demonstrate that a 
typical a-particle condensate is the Hoyle state (E x = 7.65 MeV, 0^ state in 12 C), 
which plays a crucial role for the synthesis of 12 C in the universe. The influence 
of antisymmentrization in the Hoyle state on the bosonic character of the a parti- 
cle is discussed in detail. It is shown to be weak. The bosonic aspects in the Hoyle 
state, therefore, are predominant. It is conjectured that a-particle condensate states 
also exist in heavier na nuclei, like le O, 20 Ne, etc. For instance the 0j state of lfi O 
at E x = 15.1 MeV is identified from a theoretical analysis as being a strong can- 
didate of a 4a condensate. The calculated small width (34 keV) of 0g\ consistent 
with data, lends credit to the existence of heavier Hoyle-analogue states. In non- 
self-conjugated nuclei such as n B and 13 C, we discuss candidates for the product 
states of clusters, composed of a's, triton's, and neutrons etc. The relationship of a- 
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particle condensation in finite nuclei to quartetting in symmetric nuclear matter is 
investigated with the help of an in-medium modified four-nucleon equation. A non- 
linear order parameter equation for quartet condensation is derived and solved for 
a particle condensation in infinite nuclear matter. The strong qualitative difference 
with the pairing case is pointed out. 



1 Introduction 

Cluster as well as mean-field pictures are crucial to understand the structure of light 
nuclei |Q~|[2). It is well known that many states in light nuclei as well as neutron rich 
nuclei and hypernuclei 2] have cluster structures. Recently, it was found that 
certain states in self-conjugate nuclei around the a-particle disintegration threshold 
can be described dominantly as product states of a particles, all in the lowest OS 
orbit. They are called "a-particle condensate states". Considerable theoretical and 
experimental work has been devoted to this since this idea was first put forward in 
2001 0. 

The ground state of 8 Be has a pronounced a-cluster structure ]6]|7). Its average 
density in the + ground state is, therefore, very low, only about a third of usual nu- 
clear saturation density. The two a particles are held together only by the Coulomb 
barrier and 8 Be is, therefore unstable but with a very long life time (10~ 17 s). No 
other atomic nucleus is known to have such a structure in its ground state. How- 
ever, it is demonstrated with a purely microscopic approach that, e.g. 12 C also has 
such a structure but as an excited state I181 4TT1 : the famous "Hoyle" state 11211131 . 
i.e. the 0^ state at 7.65 MeV fl4l . It is formed by three almost independent a parti- 
cles, only held together by the Coulomb barrier. It is located about 300 keV above 
the disintegration threshold into 3a particles and has a similar life time as 8 Be, i.e. 
also very long. A new-type of antisymmetrized a-particle product state wave func- 
tion, orTHSR a-cluster wave function proposed by Tohsaki, Horiuchi, Schuck, and 
Ropke l5l [l5HT8l describes well the structure of the Hoyle state. The THSR wave 
function is analogous to the (number-projected) BCS wave function [ 19], replacing, 
however, Cooper pairs by a particles (quartets). The 3a particles, to good approxi- 
mation, can be viewed to move in their own bosonic mean field where they occupy 
the lowest OS level. We, therefore, talk about an alpha particle condensate. A more 
accurate theory reveals that there exist residual correlations, mostly of the Pauli 
type, among the alpha particles and, in reality, their occupation of the OS level is 
reduced but still amounts to over 70 % ||20l - l22l . This number is typical for nuclear 
mean field approaches. The theory |5][l0] reproduces almost all measured data of 
the Hoyle state, as for instance the inelastic form factor from (<?, <?'), very accurately. 
It is predicted that the Hoyle state has about triple to quadruple volume compared 
with the one of the '"C ground state. Excitations of one alpha out of the condensate 
into 0D and 15 states of the mean field can be formed and the 2\ [21 23 24| and 
0^ ll25l states in 12 C are reproduced in this way (the latter, so far only tentatively). 
This triplet of states are precisely the ones which, even with the most modern no- 
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core shell model codes (26 , 27], cannot be reproduced at all. The new interpretation 
of the Hoyle state as an a condensate has stimulated a lot of theoretical and experi- 
mental works on a-particle condensation phenomena in light nuclei 1 23 -25 .28-37). 

The establishment of the novel aspects of the Hoyle state incited us to conjecture 
4a condensation in le O. The theoretical calculation l38l of the OCM (orthogo- 
nality condition model) type l39l succeeded in describing the structure of the first 
six + states up to about 16 MeV, including the ground state with its closed-shell 
structure, and showed that the 0^ state at 15.1 MeV around the 4a threshold is a 
strong candidate for a 4a-particle condensate, having a large a condensate frac- 
tion of 60 %. Similar gas-like states of a clusters have been predicted around their 
a cluster disintegration thresholds in self-conjugate A = An nuclei with the THSR 
wave function [5,40| and the Gross-Pitaevskii-equation approach BP . Besides the 
An nuclei, one can also expect cluster-gas states composed of alpha and triton clus- 
ters (including valence neutrons, etc.) around their cluster disintegration thresholds 
in A An nuclei, in which all clusters are in their respective OS orbits, similar to the 
Hoyle state with its (0S a ) 3 configuration. The states, thus, can be called "Hoyle- 
analogues" in non-self-conjugated nuclei. It is an intriguing subject to investigate 
whether or not Hoyle-analogue states exist mA^An nuclei, for example, 1 1 B, com- 
posed of 2a and a t cluster 114211441 or 13 C, composed of 3a and In RBT - PiTll . The 
2a +t (3a + n) OCM BUSH calculation indicates that the 1/2^ (1/2^) state at 
E x = 1 1 .95 (12. 14) MeV just above the 2a + 1 (3a + n) threshold is a candidate for 
the Hoyle-analogue. 

It has been pointed out that in homogeneous nuclear matter and asymmetric mat- 
ter a condensation is a possible phase 148-52 1 at low densities. Therefore, the above 
mentioned a-particle product states in finite nuclei is related to Bose-Einstein con- 
densation (BEC) of a particles in infinite matter. The infinite matter study used a 
four particle (quartet) generalization of the well known Thouless criterion for the 
onset of pairing as a function of density and temperature. The particular finding in 
the four nucleon case was that a-particle condensation can only occur at very low 
densities where the quartets do not overlap appreciably. This result is consistent with 
the structure of the Hoyle state as well as the 0^ state of 16 0, in which the average 
density is about one third or one fourth of the saturation density. It is interesting to 
note that the low density condition for quartetting was in the meanwhile confirmed 
in Ref. 11531 with a theoretical study in cold atom physics. 

At this point it may be worthwhile to remark that nuclear physics is predestinated 
for cluster physics. This stems from the fact that in nuclear physics there are four 
different fermions (proton-neutron, spin up-down), all attracting one another with 
about equal strength. Such a situation is very rare in interacting fermion systems. 
Most of the time there are only two species of fermions, as e.g. electrons, spin 
up-down. However, four different fermions are needed to form a quartet. This is 
easily understood in a mean field picture where the four nucleons can be put into 
the lowest 05 level of a harmonic potential, whereas were there only neutrons two 
of four neutrons would have to be put into the p-orbit which is energetically very 
penalizing, see Fig.Q] This is the reason why there is no bound state of four neutrons, 
while the a-particle is very strongly bound. However, recently experiments in cold 
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(a) (b) 




Fig. 1 (Color online) Sketch of (a) a-particle configuration with the two protons and two neutrons 
occupying the lowest 05 level in the mean field potential of harmonic oscillator shape, and (b) the 
energetically lowest configuration in the case of four neutrons with two neutrons in the 05 orbit 
and the other two in the OP orbit. 



atom physics try to trap more than one species of fermions |]54| which then also may 
open up interesting cluster physics in that field. 

The purpose of this lecture is to demonstrate the novel aspects of nuclear a- 
particle condensates, in particular, emphasizing the structure study of 12 C and lfi O 
with the THSR wave function and the OCM approach. 

The paper is organized as follows. In Sec. [2] we first review briefly the RGM 
framework to describe na nuclear states H55II56I , which is basic for the THSR wave 
function and OCM. Then, we formulate the THSR wave function and OCM. Before 
discussing the Hoyle state, we study the structure of 8 Be with the THSR wave func- 
tion, and discuss the difference between the THSR-type wave function and Brink- 
type wave function 1 57 58 1 in Sec. [3] The latter type of wave function is based on 
a geometrical, crystal-like viewpoint of the cluster structure. Section|4]is dedicated 
to a discussion of the structure of the Hoyle state, studying the antisymmetrization 
effect among the 3 a clusters, occupation probability and momentum distribution 
of a particles, and the de Broglie wave length, etc. Then, we discuss the Hoyle- 
analogue states in le O with the 4a OCM and THSR wave function, together with 
n B and 13 C. The Gross-Pitaevskii-equation approach is devoted to investigate a- 
particle condensation in heavier An nuclei. In Sec. |3J we focus on the a-particle 
condensation in nuclear matter and its relation with that in finite nuclei. The density 
dependence of the a condensation fraction is discussed and a 'gap' equation for the 
a particle order parameter is established and solved. The strong qualitative differ- 
ence with the pairing case is discussed. Finally, in Sec. [6] we present the summary 
and conclusions. 
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2 Formulation of alpha-condensation 

— THSR wave function and OCM approach — 

2.1 Resonating Group Method (RGM) 

The microscopic na wave function ¥5, a incorporating a-cluster substructures can 
in general be expressed in the following RGM form 



f»a = A|x(5)n^| = / d*%a(*)z{*), (1) 

^ 7a (a) = yi|n5^.-a 7 )nfe,|, (2) 

with A the antisymmetrizer of An nucleons. The intrinsic wave function of the i-th 
a cluster, (j) aj , is taken as a Gaussian (with size parameter b), 



<j) aj « exp 



£ (r a -r u ) 2 /(8b 2 ) 

l<k<l<A 



(3) 



representing the intrinsic spatial part of the (Os) 4 shell-model configuration, where 
{r, i, • • • ,r,-4} denote the coordinates of the four nucleons in the i-th cluster. The 
spin-isospin part in Eq. (O is not explicitly written out but supposed to be of scalar- 
isoscalar form. We will not mention it henceforth. The wave function % for the 
c.o.m. motion of the a's is chosen translationally invariant and depends only on the 
corresponding Jacobi coordinates £ = {£i,§2>" ' The function f^a (a) in 

Eq. (01 describes the a-cluster state located at the relative positions specified by a 
set of the Jacobi parameter coordinates a = {ai,a2, • • • ,a„_i}. 

The internal part of the Hamiltonian for the relevant A=4n nucleus is composed 
of kinetic energy with nucleon mass M, the Coulomb force (Vy), the ef- 

(2) (3) 

fective two-nucleon (VA ') and three-nucleon (V-j k ) interactions: 



An j-2 An An An 

^ = -E^-7b+L^+I^ , + E v®, (4) 

1=1 i<j i<j i<j<k 



where the c.o.m. kinetic energy of the total system Tq is subtracted. 
The Schrodinger equation for the fermionic na system is 

H% a =E% a . (5) 



Substituting the total wave function of Eq. into Eq. @, we obtain the equation 
of motion for the relative wave function %, 
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Jda'{H(a,a')-EN(a,a')}x(a')=0, or (M-EN) X =0, (6) 
where the Hamiltonian and norm kernels, H(a, a') and N(a,a'), are defined as 

Equation (|6]i is called the RGM equation |56l . One also can formulate the RGM 
framework for non An nuclei such as n B and 13 C with the microscopic 2a + 1 and 
3 a + n cluster model, respectively. 



2.2 THSR wave function 

In the THSR description 151 1171 . the relative wave function % in Eq. (Q~|) is expressed 
in the following na condensation form, 

X™ SR (B:R h R 2l ---,R n ) = f\(p (B:R l -X G ), (8) 

<Po(B:R) =exp(-2R 2 /B 2 ), (9) 

where R; = (r,j H +r, 4)/4 denotes the c.o.m. coordinate of the 2-th a parti- 
cle, Xq = (Ri H + R n )/n is the total c.o.m. coordinate of the na system and 

(po(B : R) represents a Gaussian with a large width parameter B which is of the nu- 
cleus' dimension. Usually, one uses Jacobi coordinates splitting off the total 
c.o.m. part of the wave function. Then the THSR ansatz for % in Eq. ([8]l is given by 

X™ SR (B : Ri,R 2 , • • • ,R„) = exp ( - 2 £ Mi |iJ , (10) 

with jij = ;/(; + 1). A slight generalization of Eq. ( TTOb is possible, taking into ac- 
count nuclear deformation (see Sec. |3j. With Eqs. © and ([1}, one can write the 
THSR wave function in the following na product form, 

%a -► (ri ) i,--- > r B)4 |THSR)=^[vra 1 V»i%"-V*], (H) 

where 

| THSR) = |THSR(B)) =A\B), (12) 

(r u ,---,r„, 4 |B) = WaxWa a ---Wa n , (13) 

where \j/ aj = (po(B : R; — Xg)^, and definitions of Eqs. ( fT2l and ( fT3l ) will be useful 
later. Equations (fTTTi^(fT3]) show the analogy of the THSR wave function with the 
number-projected BCS wave functions for pairing 



Nuclear Alpha-Particle Condensates 



7 



(r u ,---,r„, 2 |BCS) = ^^(ri^ri^^foi,^)--- ^HBr(r B) i,r„ j2 )], (14) 

where 0pair(rj,i,r»,2) denotes the Cooper pair wave function. 

The product of n identical OS wave functions in Eq. ( fTTT i reflects the boson con- 
densate character. This feature is realized as long as the action of the antisym- 
metrizer in Eq. (fl~|i is sufficiently weak. On the other hand, in the limit where B 
is taken to be B = b, the normalized THSR wave function is equivalent to an SU(3) 
shell model wave function with the lowest harmonic oscillator quanta 15911601 : for 
example, in the case of the 2a, 3a and 4a systems, they respectively are given by 

Km N 2a (B)^ a (B) = \(0s) 4 (0p) 4 ;(Xpi) = (4,0),/* = 0+), (15) 

B^b 

\im N 3a (B)V 3a (B) = |(0s) 4 (0 P ) 8 ;(Aa0 = (0,4),/* = 0+>, (16) 

B^b 

MmN 4a (B)% a {B) = \(0s) 4 (0p) n ;(Xn) = (0,0),/" = 0+), (17) 

B^b 

where the N na (B) are the normalization factors. The shell model wave functions in 
Eqs. ( fT5l l. (fT~6b and (fTTT i are the dominant configurations of the ground-state wave 
functions of 8 Be, 12 C and le O, respectively. In fact, the components of Eqs. ( fT6b 
and ( fTTI i in the ground states of 12 C and 16 have weights over 60 % and 90 %, 
respectively, because both of the states have shell-model-like compact structures. On 
the other hand, in the case of the ground state of 8 Be, the component of Eq. (TT3T > is 
as small as about 20 % but still the largest, and the remaining components distribute 
monotonously in a lot of higher SU(3) configurations, when one expands the wave 
function of the 8 Be ground state in terms of the SU(3) basis. This characteristic 
comes from a pronounced 2a cluster structure in the ground state (see Sec. |3). 

The wave functions of the quantum states in A = An nucleus can be expanded 
using the na THSR wave function, like 

^ = £A(« (m) )^a(B (m) ), (18) 

m 

where f„ a is the n a THSR wave function which has the form of 

% a (B {m) ) = yi[z™ SR (Z?( m) ;Ri,R 2 , • • ■ ,1^0^ • • • K]. (19) 

The discrete variational parameters fiW represent the generator coordinate of the 
Hill-Wheeler ansatz. The expansion coefficients fk(B^) and the corresponding 
eigenenergy for the fc-th eigenstate are obtained by solving the following Hill- 
Wheeler equation 



W^ a (fl( m >) H-E k % a (B^))f k (B^) = 0. (20) 



This equation has the same structure as the RGM equation in Eq. (|6]l but reducing it 
to a one parameter equation. The superposition of THSR wave functions fine tunes 
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the results but a single THSR wave function with an optimized B-value already 
yields excellent results as will be demonstrated below. 



2.3 na boson wave function and OCM 

In order to study the bosonic properties of the na system, one needs to map the 
microscopic (fermionic) na cluster model wave function % M in Eq. (Q]i onto an 
na boson wave function &„ a ' . The RGM framework given in Sec. 12. 11 is useful 
and appropriate for the mapping. Taking into account the normalization of ¥^ a , 
1 = (f^alffia) = (x(£)\N{£ ,£')\x(£'))' me na bosonic wave function is provided 
in the following form ll56l . 



where % represents the relative wave function with the set of Jacobi coordinates, 
% = {4i)§2)'"")§ n -i}> w i tn respect to the c.o.m. of a clusters. The square-root 
matrix N X I 2 (E,,E,') is related to the norm kernel of the na RGM wave function in 
Eq. (Q. It is noted that <pj£} depends only on the Jacobi coordinates |, and all of 

(B) 

the internal coordinates of na particles are integrated out in <P„ a . 

From the RGM equation ©, the equation of motion for is obtained in 

the form 



where "K denotes the Hamiltonian kernel defined in Eq. @. Then, one can interpret 
N-WjQT 1 / 2 as the nonlocal na boson Hamiltonian. In Eq. (22) care should be 
taken that before inversion all zero eigenvalues of the norm 3\f are properly elimi- 
nated. The eigenfunctions belonging to the zero eigenvalues are the so-called Pauli 
forbidden states Uf(r) which satisfy the condition N«f = A )IIj=i 0a,-} = 0. 

(B) 

The boson wave function has the following properties: 1) <P„ a ' is totally sym- 
metric for any 2a-particle exchange, 2) 4$ satisfies the equation motion and 
3) <&na is orthogonal to the Pauli forbidden states uf (r). In order to obtain the bo- 

(B) 

son wave function <P„ a , we need to solve the equation of motion of the bosons in 
Eq. ( l22b . Solving the boson equation, however, is difficult in general even for the 
3a case. Thus, it is requested to use more feasible frameworks for the study of the 
bosonic properties and the amount of a condensation for the Na system. One such 
framework is OCM (orthogonality condition model) 11391 . The OCM scheme, which 
is an approximation to RGM, is known to describe nicely the structure of low-lying 
states in light nuclei Il2l l2ni25l[38ll39ll63l466l . The essential properties of the na 
boson wave function <P„ a , as mentioned above, can be taken into account in OCM 
in a simple manner. We will demonstrate this below. 



^^)^l 2 X = jd^'N^X)x^) 



(21) 




(22) 
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In OCM, the a cluster is treated as a point-like particle. We approximate the non- 
local na boson Hamiltonian in Eq. (f22b by an effective (local) one, that is //(° CM ) 5 

^-i/z^-i/z^ocM) (23) 
ff (0CM)_£ ri _ rG+ £ vg(i,j)+ £ Vg(i,j,k), (24) 

i=l i<y=l i<j<k=\ 

where 7} denotes the kinetic energy of the i-th a cluster, and the center-of-mass ki- 
netic energy Tq is subtracted from the Hamiltonian. The effective local 2a and 3a 
potentials are presented as VJa (including the Coulomb potential) and V 3 e *, respec- 
tively. Then, the equation of the relative motion of the na particles with //( OCM ), 
called the OCM equation, is written as 



H 



(OCM) 



-E}*g*>=0, (25) 

(u F \4^ M) )=0, (26) 

where ujt denotes the Pauli-forbidden state of the na system as mentioned above. 
In the case of 2a system, the Pauli-forbidden states between the two a -particles are 
OS, OD and IS states with the total oscillator quanta Q less than 4. It is pointed out 
that the Pauli-forbidden states in the na system can be constructed from those of 
the 2a system ll67l . 

The bosonic property of the wave function <P can be taken into account by sym- 
metrizing the wave function with respect to any 2a-particle exchange, 

-.(OCM) -.(OCM)/. r, \ 

<&„a =S®na (1,2, •••,«), (27) 



where S denotes the symmetrization operator, S = (l/vnOLjt^*:, where the sum 
runs over all permutations 7 of the n a-particles. It is noted that the completely 
collapsed state of the na particles is forbidden within the present framework because 
of the Pauli-blocking effect in Eq. (1261) . 

The OCM equation ( f25b with the condition ( 1261 ) is solved with the help of the 
Gaussian expansion method (GEM) [68. 69]. Combining OCM and GEM provides 
a powerful tool to study the structure of light nuclei [21. 38l l44ll46ll as well as 
light hypernuclei [|4][70), because the Pauli-blocking effect among the clusters is 
properly taken into account and GEM covers an approximately complete model 
space l68ll69l . It is also useful to apply Kukulin's method f7T| for removing the 

Pauli-forbidden states Uf's from the wave function ^ > ia CM ' 1 . The present OCM- 
GEM framework, for example, in the case of ie O, can cover a model space large 
enough to describe the dilute a gas-like configuration, as well as a + 12 C cluster 
and shell-model-like ground state structures. 
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2.4 Single a-particle density matrix and occupation probabilities 

A literal interpretation of an a condensate in a finite system is that all the a parti- 
cles occupy the lowest OS-wave orbit of an a mean field potential. Due to residual 
interactions and the action of the Pauli principle, the occupation probabilities may 
spread out over several orbits, but a particular orbit should be occupied with a sig- 
nificant probability if a state is called a condensate. The occupation probability can 
be calculated by solving the eigenvalue problem of a single a-particle density ma- 
trix 1I|J22]|72H75]. 

The single a-particle density matrix for the na boson system can be defined 
with the use of the na Boson wave function 0„„ ) in Eq. (f2TT > mapped from the 
translationally invariant normalized microscopic na wave function in Eq. (UJ), 



P I ( n 1 t ) (qi,qi')=(^ T j P&{ti,tA (28) 

Pffil(«l,«l')= flidt i *$\ZuZ2"-An-l)*!£{ZlA2,"-,Z^l), (29) 
J i=2 

where q l = '^-^i — Ri — Xg is the 1st particle coordinate (Ri) with respect to 
the c.o.m coordinate of the system (Xg) and £ j denotes the relative coordinate be- 
tween the 1st particle and the remaining (n— 1) ones. The factor in Eq. (l28l is the 

Jacobian d^/dq^ Since the wave function is totally symmetric with re- 

spect to particle permutation, the choice of the 1st particle is arbitrary. The definition 
d29l ) is called the Jacobi-type one-particle density matrix. The diagonal density ma- 
trix p/ n ^(q,q) stands for the density distribution of a particles with respect to the 
c.o.m. coordinate of the na system. The eigenvalue problem of the density matrix 

P W 



(q,q')<P(q')=^<P(q) ! (30) 



gives the single a-particle orbit <p(q) and its occupation probability A, where q is 
measured from the c.o.m. coordinate of the system. The spectrum of eigenvalues of 
the density matrix gives information on the occupancy of the orbits of the sys- 
tem and it is obviously equal to that of p/ n ' t \- The occupation probability is labeled 
with the angular momentum L and the quantum number of a positive integer ni, like 
L„ L . In this article, for a single- a orbit with an angular momentum L, we denote the 
largest occupation probability as L\ («£ = 1), the second largest as L2 (ni = 2), the 
third largest as L3 («3 = 3), etc. Please notice that the positive number «£ is differ- 
ent from the number of nodes for the radial part of the corresponding single- a orbit 
<p(q) (for instance, in Fig. Qj] the single-a orbits labeled as L = and n^ = 1 (Si) 
have 2S and OS nodal behaviors for the ground state (0^) and the Hoyle state (0j ), 
respectively). 
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Let us remind that one should use the Jacobi coordinate system for the choice of 
the internal coordinates of the density matrix. If an internal coordinate system other 
than the Jacobi coordinate system is adopted (for example, that adopted by Pethick 
and Pitaevskii (PP) 1761 ). an unphysical result is obtained even for condensation of a 
finite number of ideal bosons in a harmonic trap, contrary to what PP expected 17711 . 
Two physically motivated criteria for the choice of the adequate coordinate system 
lead to a unique answer for the internal one-particle density matrix, i.e. the Jacobi- 
type internal density matrix, while the PP-type one-body density matrix does not 
satisfy the criteria (see Refs. 17411751 for details). 

In general even for the 3 a system, one encounters numerical difficulties to ob- 
tain the boson wave function mapped from the microscopic na wave function, 
<£>na = in Eq. (|2TT >. by solving the boson equation in Eq. d22l ). as mentioned 

above. Thus, it is hard in general to calculate the one-body density matrix for the a 
particle pP/ in Eq. d28l i. To overcome this difficulty, the following two approximate 
ways have so far been proposed to evaluate the density matrix. One is, as for the 
boson wave function, to use the na OCM wave function ( 1271 ) obtained by solving 
the OCM equations (|23) and ([25}, i.e. ~ <P^ M] J2TJES1- The application of 
this method was done for the 3 a and 4a systems, the results of which are presented 
in Sees. 14.1 1 and 14.21 The other is to make the following approximation for the bo- 
son wave function as proposed in Refs. (20]|22), <P ( n B a (§ ) ~ ft*/ \/WxWx)- This 
method was used for the 3a and 4a THSR wave functions. The results are discussed 
in Sec. 14.21 The two approximations give quantitative similar results for the occupa- 
tion probabilities and for the single- a orbits in the 3 a and 4a systems, which are 
obtained by solving the eigenvalue problem of Eq. (13 Oi l. 



3 THSR wave function vs Brink wave function for 8 Be 

Before discussing the Hoyle state, it is instructive to study 8 Be in some detail be- 
cause even this nucleus which is known to have intrinsically a two-alpha dumbbell 
structure |6j[7] can very well be described in the laboratory frame with the THSR 
wave function. Let us repeat Eq. (Q~|i for this particular case 

*ba=A\Z( T )Ma I \, (31) 

with the relative coordinate between the 2a particles, r = Ri — R2. Note that 
Eq. OUi is a fully antisymmetric and translationally invariant wave function in 
8 — 1=7 coordinates. Solving the RGM equation in Eq. (0) with a given Hamilto- 
nian, one obtains the energy E of 8 Be and %. The 2a boson wave function <P^ (r) 
representing the relative motion of the two a-particles, mapped from the corre- 
sponding fermionic 2a wave function fie is given in Eq. (1211) . 

*2«V) = / dr'N l l 2 {v,T')x{r'l (32) 



12 



T. Yamada, Y. Funaki, H. Horiuchi, G. Ropke, P. Schuck, and A. Tohsaki 



1.0 r THSR 




0.0 5.0 10.0 , 15.0 20.0 25.0 

rfm 



Fig. 2 (Color online) Comparison of THSR wave function with a single component "Brink" wave 
function with D — 3.45 fm (denoted by n = 1). The convergence rate with the superposition of 
several («) "Brink" wave functions is also shown. The line denoted by n = 30 corresponds to the 
full RGM solution. The Volkov No.l force is taken with Majorana parameter value M = 0.56. 
Figure is taken from Ref. 1171 . 



Expressions (f3TT > and d32l have been obtained with very high numerical accuracy 
since 50 years with excellent results for all low energy properties of 8 Be (6). The 
radial part of the 2a boson wave function r<p2a( r ) m the ground state (J K = + ) 
is shown in Fig. [2] denoted by n = 30. We see that there exist two nodes, an effect 
which stems from the Pauli principle. 

Here we will discuss two approximate forms for x( r ) : tne THSR wave function 
and the Brink cluster wave function ||57l . Let us start with the latter. In the Brink 
wave function, the a particles are placed at certain positions in space. In the case of 
8 Be, placing the 2a particles at the positions of D/2 and — D/2, respectively, this 
leads to 



Z Brink (r) = F=°exp - (R, - D/2- (R 2 +D/2) /b 



2 



,2 



P J -°exp 



(33) 



where P J=Q denotes the projection operator onto spin / = 0. Though this kind of 
geometrical, crystal-like viewpoint of the cluster structure works well for many 
cases, for instance, parity-violating 12 C+a, lfi O+a, and 40 Ca+a structures in lfi O, 
20 Ne and 44 Ti, respectively 0[2S|79]|, and also when additionally neutrons are in- 
volved ll80l . it is on the contrary known since several decades that this picture fails 
for the description of the famous Hoyle state, i.e. the 0^ state in 12 C (see Sec. 0). 
The ansatz of the two a particles being placed at a distance D from one another 
seems reasonable, since the Quantum Monte Carlo calculation with realistic two- 
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nucleon and three-nucleon potentials in Ref. |7| indeed indicates that the two a's 
are about 4 fm apart. Obviously, the parameter D can be varied to find the optimal 
position of the a-particles. The result of such a procedure is shown in Fig. |2]with 
the line denoted by n = 1 taking the optimal value D = 3.45 fm (b is kept fixed at its 
free space value, b = 1.36 fm). Qualitatively such a "Brink" wave function follows 
the full variational solution (line denoted by n = 30). However, in the outer part, 
for instance in the exponentially decaying tail quite strong differences appear. The 
squared overlap with the exact solution is 0.722. Of course,the Brink wave functions 
also can serve as a basis and it is interesting to study the convergence properties. We, 
therefore, write for the 8 Be wave function appearing in Eq. © 

%a =A\x{r)$M =£^ B c T k (r, J D« J fc), (34) 



.Brink / 



(35) 



where the indicate the various positions of the a-particles and /, are the expan- 
sion coefficients. The convergence of the squared overlap with the exact solution 
is studied where we take for the positions DW = 1 fm, = 2 fm, ■ ■ ■, D( n > = n 
fm. We start with n = 5. In Fig. [3] the convergence rate is shown as a function of 
n for the squared overlap and for the energy. The point of n = 1 is with the opti- 
mized single Brink wave function (D^ l) = 3.45 fm). We see that the convergence is 
not extremely fast but for n = 20 the squared overlap with the full RGM solution 
amounts to 0.9999. Also energy is converged to within 10~ 4 . In Fig.|2]we show the 
convergence of the 2a boson wave function r<p2a{f)- In the insert we see that there 
is still a slight change in the far tail going from n = 25 to n = 30. 

Let us now investigate the THSR ansatz for %{r). There it is assumed from the 
beginning that the a's are delocalised and a single Gaussian e ' ' centered at the 
origin with, however, a large width B 2 = b 2 + 2/3 2 , with j3 a variational parameter, 
is taken. Very much improved results over the single component Brink wave func- 
tion are obtained. With j5 = 3.24 fm the squared overlap becomes 97.24 %. How- 
ever, practically 100 % accuracy, compared with the exact solution, can be achieved 
starting with a slightly improved ansatz, i.e. with an axially symmetric deformed 
Gaussian which is then projected on the ground-state spin / = (projections on 
/ = 2,4 yield the rotational band of 8 Be) 18D . 



X THSR (r) -P^exp 



2 2 



b 2 + 2p 2 b 2 +2p 2 



exp(-r 2 /Bi) /.(fi?-^) 1 / 2 



Erf ' ' p it n, (36) 

ir V B±o ; ' 

with B 2 = b 2 + 2/3? and r 2 L = r 2 + r 2 , and Erf(x) the error function. The second line 
of Eq. d36l l is obtained from a simple calculation. 

Such an intrinsically deformed ansatz is, of course, physically motivated by the 
observation of the rotational spectrum of 8 Be indicating a large value of the corre- 
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Fig. 3 (Color online) Binding energy corresponding to f^ a (see Eq. i34\ ) with the supeiposition 
of n Brink wave functions and the squared overlap between the full RGM solution and f^ a . For 
n = 1, a single Brink wave function with optimized R = 3.45 fm is adopted. Figure is taken from 
Ref. ED. 



sponding moment of inertia. The minimization of the energy yields fi± = fix = fiy = 
1.78 fm and J3- = 7.85 fm. With these numbers, the squared overlap between the 
exact f*> a and f 9 ™ SR is with 0.9999 extremely precise. In Fig.[2]we also show that 
the THSR wave function agrees very well even far out in the tail with the "exact" 
solution with 30 "Brink" components. 

As seen above, the single component, two parameter THSR ansatz, Eq. (l36l ). 
for the relative wave function of two alpha's seems to grasp the physical situation 
extremely well. The most important part of this wave function is the outer one be- 
yond some 3 fm. There, the two alpha's are in an S wave of essentially Gaussian 
shape. The corresponding harmonic oscillator frequency is estimated to h(0 ~ 2 
MeV. Therefore, as long as the two alpha's do not overlap strongly, they swing in 
a very low frequency harmonic oscillator mode in a wide and delocalized fashion, 
reminiscent of a weakly bound gas like state. Inside the region r < 2-3 fm where 
the two alpha's heavily overlap, because of the strong action of the Pauli principle, 
the relative wave function has two nodes and small amplitude, as shown in Fig. [2] 
Contrary to the outer part of the wave function determined dynamically, the behav- 
ior of the relative wave function in this strongly overlapping region is determined 
kinematically, solely reflecting the r-dependence of the norm kernel in Eq. (JT). This 
is clearly seen from the fact that both THSR and Brink wave functions have very 
nearly the same behavior in this region. Thus, we found that the alpha's in 8 Be 
move practically as pure bosons in a relative OS state of very low frequency as long 
as they do not come into one another's way, that is as long as they do not overlap. 
One should stress that this picture holds after projection on good total momentum 
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Fig. 4 Experimental energy spectra of 12 C 1141 together with the calculated ones using the 3 a 
RGM (9). 



and good spin, that is in the laboratory frame. It is equally true, as already men- 
tioned, that in the intrinsic frame 8 Be can be described as a strongly deformed two 
alpha structure, see ansatz (l36l . reminiscent of a dumbbell. 



4 Alpha-gas like states in light nuclei 

4.1 n Ccase 

The a cluster nature of 12 C has been studied by many authors using various ap- 
proaches. Figure |4] shows the energy spectrum of 12 C lfT4l . The 0^ state, located 
near the 3 a breakup threshold, is called the Hoyle state |[T2][T3"1 . which plays an 
astrophysically crucial role in the synthesis of 12 C in the universe. Its small excita- 
tion energy of 7.65 MeV is very difficult to explain by the shell model, even using 
the most modern non-core shell model approach [26. 27 1. The fully microscopic 3 a 
cluster models |8 9|, however, succeeded in the 1970s in explaining the observed 
data such as the small excitation energy and the inelastic form factor of the (e,e') 
reaction etc., together with the structures of the ground-band states (0^ — 2\ — 4j~), 
2t, and negative-parity states (3^ — 1^). The cluster model studies with the 3a 
GCM (generator coordinate method) and 3 a RGM fU showed that the Hoyle 
state has a weakly interacting gas like 3a-cluster structure with a very large radius 
(about 1/3 of the ground-state density), whereas the ground state has a shell-model- 
like compact structure. 

This 3a gas-like nature of the Hoyle state is demonstrated in Fig. [5] in which the 
overlap between a Brink-type wave function and the full RGM solution obtained by 
solving the 3 a RGM equation in Eq. © is shown. The overlap is quite poor and 
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Fig. 5 Structure of the Oj 
state shown by the overlap be- 
tween the Brink-type cluster 
wave function of the isosceles 
configuration and the exact 

wave function. Figure 
adopted from Ref. t2llfl . 



Table 1 Comparison of the total energies, r.m.s. radii (/? r .m.s.)> an d monopole strengths (M(0j — ► 
0^ ) ) for 12 C given by solving Hill-Wheeler equation based on Eq. ifJJ and by Ref. (9) . The effective 
two-nucleon force Volkov No. 2 182 1 was adopted in the two cases for which the 3 a threshold 
energy is calculated to be —82.04 MeV. 







THSR w.f. 










(Hill-Wheeler) 


3a RGM (9) 


Exp. 






-89.52 


-89.4 


-92.2 


E (MeV) 


0+ 


-81.79 


-81.7 


-84.6 


Sr.m.s. (fm) 


ot 


2.40 


2.40 


2.44 




3.83 


3.47 




M{0+ ->0j") (fm 2 ) 




6.45 


6.7 


5.4 



in the best case the squared overlap reaches only about 50%. This means that the 
0^ state has a distinct clustering and has no definite spacial or geometrical config- 
uration. The situation is also pointed out in a recent work itTTl . in which about 55 
components of the Brink-type wave functions are needed to reproduce accurately 
the full RGM solution for the Hoyle state. However, this Hoyle-state wave function 
is shown to be almost completely equivalent to a "single THSR wave function" as 
discussed in next section. 



4.1.1 THSR description of the Hoyle state 

The total wave function for 12 C in the THSR description is obtained by solving the 
Hill-Wheeler equation based on Eqs. (U~8T i. (TT9b and ( f20b . TableQ]shows the results of 



Nuclear Alpha-Particle Condensates 



17 



Tl 




10' 







0.5 



1.5 



2 



2.5 



3 



3.5 



q [fm 1 ] 



Fig. 6 (Color online) Comparison of the experimental inelastic form factor of l2 C(e,e') with the 
RGM (denoted by cluster), THSR (BEC) and FMD calculations. Figure is adopted from Ref. fill . 



the energies, r.m.s. radii, and monopole strengths in the THSR description together 
with those of the full 3 a RGM calculation and the data. One can see that the THSR 
description succeeds to reproduce the properties of the two + states. Inspecting the 
r.m.s. radii, the Hoyle state has a volume 3 to 4 times larger than that of the ground 
state of 12 C. The inelastic form factor of 12 C from the ground state to the Hoyle 
state in the THSR description is displayed in Fig. [6] We reproduce very accurately 
the experimental data. 

In order to study how good a single 3 a THSR wave function reproduces the full 
RGM solutions, we use the THSR wave function with axially symmetric deforma- 
tion, presented as 



where | j 2 are the two Jacobi coordinates with fi\ = 1/2 and jj.2 = 2/3, and fi± and 
J3 Z are the deformation parameters with f3± = fix = fiy- The wave function with good 
total spin J — is written as 



f-5a(Ai,/y =A{ * 3 ™ SR ( J3 ± , ft ) 0„ 0„ a } , 



(37) 




(38) 



(39) 
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Fig. 7 (Color online) Contour map of the energy surface in the two parameter space (P±,p z ) for 
(left) <Pf= (/3j_,/3,) in Eq. ([39) and for (right) 1 ff= (/3 ± , J 8 ; ) in Eq. go) orthogonal to the ground 
state. 



where P J is the angular momentum projection operator. In what concerns the THSR 
wave function for the description of the Hoyle state, the situation is slightly more 
complicated than in the 8 Be case by the fact that the loosely bound 3a configuration 
is now no longer the ground state but the Oj state at 7.65 MeV excitation energy (as 
a side remark, let us mention that usually Bose-Einstein condensates of cold atoms 
also are not the ground states of the systems which are given by small crystals) . As 
discussed in Sec. the wave function (T39t has the dominant configuration of the 
ground state in the limit of j5± = j3- = 0. Thus, in order to discuss the Hoyle state, 
we have to use the 3 a wave function which is orthogonal to the ground state, 
expressed as 

^°05 ±J ft) =P J=0 Pl s V3 a (P±,p z ), (40) 

where P|' s keeps the wave function in Eq. d40b to be orthogonal to the ground-state 
wave function, i.e. P|' s = 1 — |0j~)(0^~|. 

On the left side of Fig. [7] we show the contour map of the energy surface corre- 
sponding to the state (f39t in the two parameter space (j8j_, j3 z ), defined as 

F ,a a , , 4n 

WuW (^(/L,ft)|s£O0 x> ft)> ' ( } 

where H is the microscopic Hamiltonian of 12 C used in the 3 a RGM calculation. 
One sees a minimum at jSj_ = 1.5 fm and J3 Z = 1.5 fm, which means a spherical 
shape. The minimum energy of —87.68 MeV is about 1.7 MeV higher than the 
total energy of -89.4 MeV obtained by the full 3a RGM calculation (see Table[T])- 
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Fig. 8 (Color online) Contour map of the squared overlap of the normalized THSR wave function 
< f^°(j8 ± ,j8j) in Eq. (|40j, orthogonal to the ground state, with the full RGM solution. 



When the Hill-Wheeler equation in Eq. (|20T > is solved in the two-parameter space of 
f5± and /3 Z , we can reproduce the total energy of the RGM result. 

On the right side of Fig. [7] the contour map of the energy surface correspond- 
ing to the state d40b orthogonal to the ground state is displayed, where we use the 
ground-state solution of the Hill-Wheeler equation in the two-parameter space of 
j3x and j5 z . We see an energy minimum at j5± = 5.2 fm and fi z = 1 .5 fm in the pro- 
late region of the map and a second energy minimum at j5± = 2.6 fm and j3- = 7.5 
fm in the oblate region. The minimum energy value is —81.75 MeV. This value is 
almost the same as the total energy of —81.67 MeV obtained by the full 3a RGM 
(see Table[T]i- The minimum energy of —8 1 .75 MeV is close to the second minimum 
energy of —81.67 MeV, and there is a valley with an almost fiat bottom connecting 
these two minima. This means that the energy of the spherical configuration is only 
slightly higher than that of the deformed configuration, that is, the energy gain due 
to the deformation is small. 

A very remarkable result from the right side of Fig. [7]is that the wave function 
at the minimum energy point (j3x = 5.3 fm and j3 z = 1.5 fm) has 99.3 % squared 
overlap with the full RGM solution (see Fig. HJ, although the spherical wave func- 
tion (/3x = Pz = 4-0 fm) gives already a squared overlap of 92 %. The THSR wave 
function Eq. d37l i is of Gaussian type with a wide extension, centered at the origin. 
It is completely different from a Brink type wave function with the three a -particles 
placed at definite values in space. A slight improvement of Eq. d40b can still be 
achieved in taking the j3, parameters as Hill- Wheeler coordinates and superpose a 
couple of wave functions of the type d40b with different width parameters. Practi- 
cally 100% squared overlap with the wave function of the full RGM result is then 
achieved. It should be pointed out that the superposition of several Gaussians of the 
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Fig. 9 (Color online) Pictorial representation of the THSR wave function for n = 3 ( 12 C). The three 
a-particles are trapped in the OS-state of a wide harmonic oscillator (B) and the four nucleons of 
each a are confined in the 0.s-state of a narrow one (b). All nucleons are antisymmetrized. 



type (|40b does not at all change the physical content of the THSR wave function as 
a wide extended distribution centered around the origin. Therefore, the Hoyle state 
can be seen as three almost inert a-particles moving in their own mean field poten- 
tial, to good approximation given by a wide harmonic oscillator, whereas the a's are 
represented by four nucleons captured in narrow harmonic potentials. The situation 
is given as a cartoon in Fig. [9] 



4.1.2 Influence of antisymmetrization and orthogonalization 

A crucial question is whether for the Hoyle state the THSR wave function ([TJ with 
^ can be considered to good approximation as a product state of a particles con- 
densed with their c.o.m. motion into the 05 orbital. For this, one has to quantify the 
influence of the antisymmetrizer A in Eq. (fTJ. A direct way to measure the influence 
of antisymmetrization is to study the following expectation value of the antisym- 
metrizer A, 

(B\A\B) 

where \B) is the THSR wave function in Eq. ( PT2b without the antisymmetrization, 
that is, just the product state y/ ai \lfa 2 Va 3 in Eq. H3i . The normalization of the an- 
tisymmetrizer A is chosen so that N(B) becomes unity in the limit where the inter- 
cluster overlap disappears, i.e. for the with parameter B — > °°. 

The result of N(B) is shown in Fig. [10] as a function of the width parameter 
B. We chose, as optimal values of B for describing the ground and Hoyle states, 
B = B g = 2.5 fm and B = Bh = 6.8 fm, for which the normalized THSR wave 
functions give the best approximation of the ground state 0^ and the Hoyle state 
Oj, respectively, which are obtained by solving the Hill-Wheeler equation (l20l . 
The squared overlaps are 0.93 and 0.78, respectively. From Fig. [l0]we find that 
N(B H ) — 0.62 and N{B g ) ~ 0.007. These results indicate that the influence of the 
antisymmetrization is strongly reduced in the Hoyle state compared with the in- 
fluence in the ground state. An important point in the present consideration is that 
the THSR wave function at B = Bh is not automatically orthogonal to the ground 
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Fig. 10 (Color online) Expectation value of the antisymmetrization operator for the product state 
|B). The value at the optimal B values, B g for the ground state and Bh for the Hoyle state, are 
denoted by a circle and a cross, respectively. 
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Fig. 11 (Color online) Energy curve in the orthogonal space to the ground state, denoted by Ep(B), 
together with E(B). The values at the optimal B values, B s and B# for the ground state and Hoyle 
state, respectively, are marked by a circle and a cross, respectively. 
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state. This is contrary to the situation with condensed cold bosonic atoms, for which 
the density is so low that the overlap of the electron clouds can, on average, be to- 
tally neglected. In the present case, the squared overlap of |THSR(Z? = Bh)) with 
|THSR(B = B g )) (or with the ground state 0|" obtained by solving the Hill-Wheeler 
equation) is less than 0.12. This small value indicates that the orthogonality with the 
ground state is nearly realized. 

An explicit orthogonalization with |THSR(B)) to the ground state 0^ obtained 
by solving the Hill- Wheeler equation gives non-negligible effects for a quantita- 
tive description of the Hoyle state with the THSR wave function. As mentioned in 
Sec. 14.1.11 the normalized THSR wave function orthogonal to the ground state 0^ 
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(~ p| g s) |THSR(B))) gives a squared overlap of 0.92 (for B = 6.1 fm) with the 
state obtained by solving the Hill-Wheeler equation, although the squared overlap 
using the normalized THSR wave function without the orthogonalization gives al- 
ready a value of 0.78. In addition, as shown in Fig. [TT] the energy curves for the 
THSR wave function, 

(THSR(fl)|//|THSR(fl)) 
1 ' (THSR(B)|THSR(B)) ' ( ' 

indicates a minimum corresponding to the ground state at B ~ B g , but the second 
minimum corresponding to the Hoyle state is not present. This is due to the fact 
that the THSR state with B = B H , |THSR(B = B H )), still includes the ground-state 
component of about 10 %, as mentioned above. In fact, if one calculates the energy 
taking into account the explicit orthogonalization to the ground state, 

= (p{ g - s) THSR( j B)|//|p| g - s) THSR(fi)) 
P (p| s ' s) THSR(B)|p} g - s) THSR(B)) 

there appears the minimum corresponding to the Hoyle state at B ~ Bh, as shown 
in Fig.QT] Thus, the small admixture of the ground-state components to the Hoyle 
state is never negligible, and explicit elimination by P]° plays an essential role 
to describe the Hoyle state. It is true that the effect of the antisymmetrization is 
not negligible even for the Hoyle state in the sense that the projection operator P^' 
excludes the compact ground-state components which are strongly subject to the an- 
tisymmetrizer. Nevertheless, it is worth emphasizing that as a result of the explicit 
orthogonalization to the ground state, the Hoyle state can not have a compact struc- 
ture but has a dilute density, for which, in the end, the effect of antisymmetrization 
is small. 



4.1.3 Alpha-particle occupation probabilities, momentum distribution, and 
the de Broglie wave length in the Hoyle state 

Direct quantities indicating how well the Hoyle state is described by a product state 
of three a's are the a-particle occupation probabilities and single particle orbits, 
which are obtained by diagonalizing the internal single a-particle density matrix 
Pint diiQ') defined in Eq. ( f28l . The occupation of the single-a orbits of the Hoyle 
state is shown in Fig. [12] One finds that the a particles occupy the S\ orbit to over 
70 %, and those for other orbits are very small. This means that each of the three 
a particles in the 0^ state is in the S\ orbit with occupation probability as large as 
about 70 %. The radial behavior of the S\ orbit is illustrated with the solid line in 
Fig.QjJb). We see no nodal behavior but small oscillations in the inner region (r < 4 
fm) and a long tail up to r ~ 10 fm. For reference, the radial behavior of the 5- wave 
Gaussian function, <po s (r) = Nq s (B) exp(— r 2 / (2B 2 )), is drawn with the dashed line 
in Fig. [T3l b). where the size parameter B is chosen to be 3.6 fm, and Nq s (B) denotes 
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Fig. 12 (Color online) Occupation of the single-a orbitals of the Hoyle state of 12 C compared with 
the ground state 



the normalization factor. The radial behavior of the S\ orbit is similar to that of 
the 5-wave Gaussian function, in particular, in the outer region (r > 4 fm), whereas 
a slight oscillation of the former around the latter can be seen in the inner region 
(r < 4 fm). Thus, the Hoyle state can be described as the product state of (OS)^ 
being realized with a probability of over 70 %. 

In the case of the ground state of 12 C, the a-particle occupations are equally 
shared between Si, D\ and G\ orbits (see Fig. [T2l. thus invalidating a condensate 
picture for the ground state. These occupancies can be explained quite well from the 
following fact: The ground state has as main configuration the SU(3) shell model 
wave function = (04). Figure \l3\ a) demonstrates the radial parts for the Si-, 

D\- and Gi-orbits, the numberof nodes of which are two, one and zero, respectively. 
Reflecting the SU(3) character, the radial behavior of the three orbits is similar to 
those of the harmonic oscillator wave functions {unl) with Q = 4, U02, M21 and M40, 
respectively, where (L) denotes the number of nodes (orbital angular momentum). 
We see that the radial parts of the single a-particle orbits oscillate strongly in the 
inside region (r < 4 fm). This is due to the important Pauli blocking effect for the 
ground state with its compact shell-model-like structure. 

Another important quantity to demonstrate the 3 a condensate nature of the Hoyle 
state is the momentum distribution of a single-a particle. It is defined as a double 
Fourier transformation of the internal single-a density matrix p^^,^') defined in 
Eq. (E8]l, 



Let us remind that p (k) would have a 5-function like peak around k = for an ideal 
dilute condensed state in homogeneous infinite matter. 




(45) 
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Fig. 13 Radial parts of the single a orbits, (a) Si (solid line), D\ (dashed) and G\ (dotted), in the 
0^ state, and (b) the Si (solid) orbit in the 0^ state compared with an S-wave Gaussian function 
(dotted), rtpoj, with the size parameter B = 3.6 fm (see text) 1211 . Note that all the radial parts in 
figures are multiplied by r. 



The momentum distributions of the a particle, p(k) and k 2 x p(k), are shown 
for the 0^ and 0^ states in Fig. [14] Reflecting the dilute structure of the Hoyle, 
we see a strong concentration of the momentum distribution in the k < 1 fm 
region, and the behavior of p (k) is of the (5-function type, similar to the momentum 
distribution of the dilute neutral atomic condensate states at very low temperature 
trapped by an external magnetic field lf83l . On the other hand, the ground state 
has higher momentum components up to k ~ 6 fm~ 1 as seen from the behavior of 
k 2 x p (k) reflecting the compact structure. The above results for the radial behavior 
of the S\ orbit, occupation probability and momentum distribution for the state 
again lead us to conclude that this state is of the 3 a condensate character with as 
much as about 70% occupation probability. 

The de Broglie wave length of the a's moving in the Hoyle state is an interest- 
ing quantity. It can be estimated from the resonance energy of 8 Be being roughly 
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Fig. 14 Momentum distribution of the a particle in 12 C, (a) p(k) and (b) k 2 p(k), for the 0^ (solid 
line) and 0j (dotted) states HTl . 



100 keV. Otherwise, one can estimate the kinetic energy of the a-particles from a 
bosonic mean field picture using the Gross-Pitaevskii equation |j4~T| (see Sec. 14.31 ). 
The mean field potential of a-particles in the Hoyle state (see Fig. l23l indicates 
the position of the single a particle energy (180 keV). The kinetic energy of the 
single a particle is calculated to be 380 keV. From this, the de Broglie wave length 
A = 2n{=yr L E a )~' i l 2 is, therefore, estimated to be of a lower limit of approximately 
20 fm. A more reliable estimate of the de Broglie wave length is to use the expecta- 
tion value of k 2 for the wave number k of the a particle in the Hoyle state, evaluated 
from the momentum distribution of the alpha particle, p (k), in Fig. [14] obtained by a 
3a OCM calculation ED- The result is A = 2n/^J{k 2 ) — 20 fm, consistent with the 
previous value. These estimates all indicate that the de Broglie wave length is much 
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Fig. 15 (Color online) Theoretical interpretation of the 0^,2^ and Oj states. 



longer than the inter a-particle distance, contrary to what is claimed in Ref. ||84l . 
using qualitative arguments. 

4.1.4 Family of the Hoyle state: 2+ and 0^ 

In the previous section, we found that the Hoyle state has a dilute 3a-condensate- 
like structure with a main configuration of (OS) ^. Then, an excited state of the Hoyle 
state, for example, a 2 + state with (0S)^(0D) a , may exist somewhat higher up in 
energy than the Hoyle state. Itoh et al. observed the 2t state at 2.6 ± 0.3 MeV above 
the 3a threshold with a width of 1 .0 ± 0.3 MeV by measuring a particles decaying 
from excited 12 C states with inelastic a scattering (28]. This state was quite recently 
confirmed by experiment with a high-energy-resolution magnetic spectrometer ||30l . 

A deformed calculation using the THSR wave function in Eqs. d37l i and ( f38T > was 
performed for the 2 + state of 12 C. Projecting on good angular momentum with a 
treatment of resonances yields the position of the 2j-state in 12 C (2.1 MeV above 
3a threshold) which is in good agreement with the experimental value |23l 28l. Also 
the calculated width (0.64 MeV) gives a quite reasonable estimate of the data. De- 
tailed investigation of the wave function of the 2^ -state shows that it can essentially 
be described in lifting out of the condensate state with the three a's in the OS- 
orbit, one a-particle in the next OD-orbit. It is tempting to imagine that the 0^-state 
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which, experimentally, is almost degenerate with the 2^ -state, is obtained by lifting 
one a-particle into the lS-orbit. Preliminary theoretical studies ||251 indicate that 
this scenario might indeed apply. However, the width of the 0^ state is very broad 
(~ 3 MeV), rendering a theoretical treatment rather delicate. Further investigations 
are necessary to validate or reject this picture which is shown graphically in Fig.lTSl 
Anyway, it would be quite satisfying, if the triplet of states, (0 \ , 2%, 0t) could all be 
explained from the a-particle perspective, since those three states are precisely the 
ones which cannot be explained within a (no core) shell model approach H26II27I . 

4.1.5 Precursors of a 3a condensate state: 3j~ and lj" 

The 3~ state at E^ a = 2.37 MeV measured from the 3a threshold in Fig. [4]is inter- 
esting from the point of view of the dilute a condensation. If the state is a condensate 
with all of the 3 a particles in the P orbit, there is the possibility of a superfuid with 
a vortex line, similar to the rotating dilute atomic condensate at very low temper- 
ature [83). Thus, it is an intriguing problem to study the structure within the 3 a 
OCM. The OCM lOTI reproduces well the energy of the 3j~ state as well as the lj~ 
state (Ej, a = 3.57 MeV) with respect to the 3a threshold. 

The calculated nuclear radius of the 3~ state is 2.95 fm, the value of which is 
larger than that for the ground state (0 while it is smaller than that for the J state. 
This suggests that the structure of the 3~ state is intermediate between the shell- 
model-like compact structure (0^) and the dilute 3a structure (0^). The occupation 
probabilities of the single-a orbits of the state are 44.7% for Pi -orbit and 27.9% for 
F\ -orbit. Although the concentration of the single Pi orbit amounts to about 50%, the 
radial behavior of the single-a orbit has two nodes in the inner region. However, the 
amplitude of the inner oscillations is significantly smaller than that for the ground 
state in Fig.[T3la) 12TI . The small oscillations indicate a weak Pauli-blocking effect, 
and thus, we can see the precursor of a 3a condensate state ED , although the 3~ 
state is not an ideal rotating dilute 3 a condensate. 

As for the lj~ state, the calculated nuclear radius, 3.32 fm, is larger than that of 
the ground state and the 3]~ state but is still smaller than that of the 0^ one. The occu- 
pation probabilities of the a particles in the 1 J" state are 35 % for Pi orbit and 16 % 
for Pi orbit. Thus, there is no concentration of the occupation probability to a single 
orbit like in the 0^ state. Since the a particles in the state are distributed over 
several orbits, the state is not of the dilute a-condensate type. On the other hand, 
the radial behavior of the Pi orbit has two nodes in the inner region, the behavior 
of which is rather similar to the 2P harmonic oscillator wave function. However, 
the Pi orbit has a P-wave Gaussian-type behavior. Also the oscillatory behavior of 
the Pi orbit for < r < 2 fm is similar to the one of the Si orbit in the 0^ state in 
Fig. Qj] These interesting behaviors of the Pi orbit indicate some signal of dilute a 
condensation, reflecting the relatively large nuclear radius (3.32 fm) of the 1 j~ state. 
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In the previous section, we showed that the Hoyle state, which has about one third 
of saturation density, can be described, to good approximation, as a product state 
of three a-particles, condensed, with their c.o.m. motion, into the lowest a mean 
field OS-orbit |2() 21 1. These novel aspects indicate that the Hoyle state has a 3 a- 
condensate-like structure. The establishment of the novel aspect of the Hoyle state 
naturally leads us to speculate about the possibility of 4a-particle condensation in 

16 o. 

The experimental + spectrum of lfi O up to about the 4a disintegration threshold 
is shown in Fig. [16] In the past, the 0+ (g.s), 0^ and 0^ states up to about 13 MeV 
excitation energy has very well been reproduced with a semi-microscopic cluster 
model, i.e. the a+ 12 C OCM (Orthogonality Condition Model) J64)- In particular, 
this model calculation, as well as that of an a+ 12 C GCM (Generator-Coordinate- 
Method) one |85), demonstrates that the 0^ state at 6.05 MeV and the 0^ state 
at 12.05 MeV have a+ 12 C structures [86] where the a-particle orbits around the 
12 C(0 1 h )-core in an 5-wave and around the 12 C(2 1 h )-core in a D-wave, respectively. 
Consistent results were later obtained by the 4a OCM calculation within the har- 
monic oscillator basis J66] . However, the model space adopted in Refs. [64 66 85 1 is 
not sufficient to account simultaneously for the a+ 12 C and the 4a gas-like config- 
urations. On the other hand, the 4a-particle condensate state was first investigated 
in Ref. (5) and its existence was predicted around the 4a threshold with the THSR 
wave function. While the THSR wave function can well describe the dilute a clus- 
ter states as well as shell model like ground states, other structures such as a+ 12 C 
clustering can not be treated and are only incorporated in an average way. Thus, it 
is important to explore the 4a condensate without any a priori assumption with re- 
spect to the structure of the 4a system. For this purpose, a full four-body 4a OCM 
calculation with Gaussian basis functions was performed ||3~8"1 . This model space is 
large enough to cover the 4a gas, the a+ 12 C cluster, as well as the shell-model con- 
figurations. In this section, we first present the results of the 4 a OCM calculation, 
and then discuss a recent analysis with the THSR wave function for the 4a system. 

4.2.1 4a OCM analysis 

The 4a OCM Hamiltonian was presented in Eq. d24b . The effective a-a interaction 
is constructed by the folding procedure from an effective two-nucleon force ll87l 
including the Coulomb interaction. One should note that the folded a-a potential 
reproduces the a-a scattering phase shifts and energies of the 8 Be ground state and 
of the Hoyle state. The three-body force V^jJ was phenomenologically introduced so 
as to fit the ground state energy of 12 C. In addition, the phenomenological four-body 
force Via was adjusted to the ground state energy of 16 0. The origin of the three- 
body and four-body forces is considered to be deducible from the state dependence 
of the effective nucleon-nucleon interaction and the additional Pauli repulsion be- 
tween more than two a-particles. However, they are short-range, and hence only act 
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Fig. 16 Comparison of en- 
ergy spectra among exper- 
iment, the 4a OCM calcu- 
lation (38), and the THSR 
treatment (22|. Dotted line 
denotes the 4a threshold. Ex- 
perimental data are taken from 
Ref . {U\ and from Ref . [35) 
for the Of state. 
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Table 2 Energies (E — £4^), r.m.s. radii (R), and monopole transition matrix elements to the 
ground state [M(E0)]. S exp . and M(E0) exp are the corresponding experimental data. The finite 
size of a particle is taken into account in R and M(E0) (see Ref. (21 1 for details). 
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in compact configurations. The expectation values of those forces is less than 10 % 
of the one of the corresponding two-body term. 

The J K = + energy spectrum obtained by the A a OCM is shown in Fig. [16] 
We can reproduce the full spectrum of + states up to about the 4a disintegration 
threshold, and tentatively make a one-to-one correspondence of those states with the 
six lowest + states of the experimental spectrum. In view of the complexity of the 
situation, the agreement is considered to be very satisfactory. 

We show in Table [2] the calculated r.m.s. radii and monopole transition matrix 
elements to the ground state, together with the corresponding experimental values. 
The r.m.s. radius of the ground state is reproduced well, and those for the other five 
+ states are by about 10 % or more larger than the ground state. The M(E0) values 
for the 0^ and 0^ states are consistent with the corresponding experimental values. 
The M(E0) value for the 0^ state is accurate only within a factor of two. 
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Fig. 17 Spectroscopic factors of the a+ 12 C(L£) channels (Lj = 0|, lj ,2|,3j ,4^,0^) in the six 
0+ sates of 16 0. 



In order to analyze the obtained wave functions, it is useful to study the overlap 
amplitude y(r) and spectroscopic factor S 2 , which are defined as follows: 



S(r'-r) 



y L (r')^ L ( 12 C) 



|f(0 6 + ) 



S z = / dr[ry(r)]\ 
Jo 



(46) 
(47) 



where <P L ( l2 C) is the wave function of 12 C, given by the 3a OCM calculation [21 ], 
and r is the relative distance between the center-of-mass of 12 C and the a parti- 
cle. From this quantity we can see how large is the component in a certain a+ 12 C 
channel which is contained in the wave functions obtained by the 4a OCM. The 
results of S 2 factors are shown in Fig. [T7] Since the ground state has a closed shell 
structure with the dominant component of SU(3)(A,ju) = (0,0), the values of the S 2 
factors for 0^ in Fig.[l7]can be explained by the SU(3) nature of the state. As men- 
tioned above, the structures of the 0^ and Ot states are well established as having 
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Fig. 18 (Colors online) Overlap amplitudes multiplied by r, defined by Eq. )46| l, for the 0^ state 
in 16 0. 



O£+ 12 C(0j~) and a+ i2 C(2~\) cluster structures, respectively. These structures of the 
0^ and 0^ states are confirmed by the 4a OCM calculation. In fact, one sees that 
the S 2 factors for the a+ 12 C(0|) and a+ 12 C(2 1 h ) channels are dominant in the 
and 0^ states, respectively. 

On the contrary, the structures of the observed O4 , 0^ and 0g states in Fig. [16] 
have, in the past, not clearly been understood, since they have never been discussed 
with the previous cluster model calculations [64. 66. 85 1. Although Ref. Q predicts 
the 4a condensate state around the 4a threshold, it is not clear to which of those 
states it corresponds to. As shown in Fig. Q~6] the 4a OCM calculation succeeded, 
for the first time, to reproduce the O4 , 0^ and 0g~ states, together with the 0^, 0^ 
and Ot states. From the analyses of the overlap amplitudes and the S 2 factors (see 
Fig. [nil, the 4a OCM showed that the 0| and 0^ states mainly have a+ 12 C(0^) 
structure with higher nodal behavior and an a+ 12 C(l~) structure, respectively. The 
monopole strength of the 0^ state is reproduced nicely within the experimental error. 

In Table[2] the largest r.m.s. radius is about 5 fm for the 0^ state. Comparing with 
the relatively smaller r.m.s. radii of the 0^ and 0^ states, this large size suggests 
that the 0^ state may be composed of a weakly interacting gas of a particles of 
the condensate type. In addition, the 0g state has a large overlap amplitude with 
the a+ 12 C(0^) channel with a S 2 factor of about 0.5 (see Fig. [TTt . whereas the 
amplitudes in the other channels are much suppressed (see Fig.[T8l. The amplitude 
in the Hoyle-state channel has no oscillations and a long tail stretches out to ~ 20 
fm. 

While a large size is generally necessary for forming an a condensate, the best 
way for its identification is to investigate the single- a orbit and its occupation prob- 
ability, which can be obtained by diagonalizing the one-body (a) density matrix as 
defined in Eq. ([28]) 120112 Ill72ti75ll . As a result of the calculation of the L = case, a 
large occupation probability of 61 % of the lowest OS-orbit is found for the Og^ state, 
whereas the other five + states all have appreciably smaller values, at most 25 % 
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(Oj )■ The corresponding single-a S orbit is shown in Fig. [19] It has a strong spatially 
extended behavior without any node (OS). This behavior is very similar to that of the 
overlap amplitude of Og for the a+ 12 C(0j) channel shown in Fig. [18] These results 
indicate that a particles are condensed into a very dilute OS single-a orbit, see also 
Refs. 12111881 . In addition, Figure l20l shows the momentum distribution k 2 p(k) of 
a particles in the six + states defined in Eq. (|43T >. One sees that the momentum 
distribution of the 0^ state concentrates strongly in the narrow region k < 1 fm _1 , 
and the behavior is quite similar to that of the Hoyle state in Fig. [14] Thus, the Og 
state clearly has, according to our calculation, a 4a condensate character. 

Comparing the single-a orbit of the Og state in Fig. [19] with that of the Hoyle 
state shown in Fig. [13] one can see an almost identical shape. This is also an im- 
portant indication that the Og state has a-particle condensate nature. Of course, the 
extension is slightly different because of the smallness of the system. The nodeless 
character of the wave function is very pronounced and only some oscillations with 
small amplitude are present in 12 C, reflecting the weak influence of the Pauli prin- 
ciple between the a's, as discussed in Sec. 14. II On the contrary, due to the much 
reduced ground-state radii, the "a-like" clusters strongly overlap in 12 C and ie O, 
producing strong amplitude oscillations which take care of antisymmetrization be- 
tween clusters 12TII381 . In fact, on sees in Fig. [19] that the single-a orbit for the 
ground state has maximum amplitude at around 3 fm and oscillations in the interior 
with two nodal (2S) behavior, due to the Pauli principle and reflecting the shell- 
model configuration (also see Fig.[l3]for the ground state of 12 C). 

The a decay width constitutes a very important information to identify the Og 
state from the experimental point of view. The width II can be estimated, based on 
the /^-matrix theory 11891 , 
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0+ states of ls O. 



Table 3 Partial a widths in 16 0* decaying into possible channels and the total width. The reduced 
widths defined in Eq. d46t are also shown, a is the channel radius. 



l2 C{Q+) + a l2 C(2+) + a n C(0+) + a 

(a = 8.0 fm) (a = 7.4 fm) (a = 8.0 fm) Total 



r L (keV) 



26 
0.006 



0.004 



2x 10" 
0.15 



34 



r L = 2P L {a)-Yt(a), 

L{a) = Fl(ka) + G 2 L {kaY 
^(a) = el(a)y^(a), 

^( fl ) = i^, 9l(a) = ^l(a), (48) 

where k, a and ji are the wave number of the relative motion, the channel radius, 
and the reduced mass, respectively, and Fl, Gl, and Pl{o) are the regular and irregu- 
lar Coulomb wave functions and the corresponding penetration factor, respectively. 
The reduced width of 9% (a) is related with the overlap amplitude defined in 
Eq. i46i . In Table[3] we show the partial a decay widths II of the 0^ state decaying 
into the a-r- 12 C(0^), a+ 12 C(2j) and a+ 12 C(0j) channels, and also the total a 
decay width which is obtained as a sum of the partial widths, and reduced widths 
0l(a) defined in Eq. (|46| |. Experimental values are all taken as given by the decay 
energies. Thus the excitation energy of the calculated 0^ state is assumed to be the 
experimental value, i.e. 15.1 MeV. 

The obtained very small total a decay width of 34 keV, in reasonable agreement 
with the corresponding experimental value of 160 keV, indicates that this state is 
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unusually long lived. The reason of this fact can be explained in terms of the present 
analysis as follows: Since this state has a very exotic structure composed of gas-like 
four alpha particles, the overlap between this state and a + 12 (1(0^) or a + 12 C(2^) 
wave functions with a certain channel radius becomes very small, as this is, indeed, 
indicated by small Ol(a) values, 0.006 and 0.004, respectively, and therefore by 
small "/1(a) values. These largely suppress the decay widths expressed by Eq. ( |48T > 
in spite of the large values of the penetration factors caused by large decay ener- 
gies 7.9 MeV and 3.5 MeV into these two channels, a + n C(0j) and a + 12 C(2^), 
respectively. On the other hand, the decay into the a + 12 C(0j) channel is also sup- 
pressed due to the very small penetration caused by the very small decay energy 
0.28 MeV, even though the corresponding reduced width takes a relatively large 
value 9l(a) = 0.15. This is natural since the 0^ state of 12 C has a gas-like three- 
alpha-particle structure. It is very likely that the above mechanism holds generally 
for the alpha condensate states in heavier na systems, and therefore the alpha con- 
densate states can also be expected to exist in heavier systems as a relatively long 
lived resonance. 

As for the decay widths of the O4 and 0^ states, as evaluated by Eq. ( |48] i, they are 
shown in Table|2] The calculated width of the 0^ state is ~ 150 keV, which is quite a 
bit larger than that found for the 0^ state ~ 50 keV. Both are qualitatively consistent 
with the corresponding experimental data, 600 keV and 185 keV, respectively. We 
should note that our calculation consistently reproduces the ratio of the widths of the 
O4 , 0^, and 0g states, i.e. about 3:1:1, respectively, though the magnitudes of the 
widths are underestimated by about a factor of four with respect to the experimental 
values (see Table |2J. The reason why the width of the O4 state is larger than that 
of the 0^ state, though the 0^ state has lower excitation energy, is due to the fact 
that the former has a much larger component of the a+ 12 C(0 1 F ) decay channel, 
reflecting the characteristic structure of the 0j~ state. The 4a condensate state, thus, 
should not be assigned to the 0^ or 0^ state 12211381 but very likely to the 0g state. 

4.2.2 THSR wave function analysis 

As mentioned already, the first investigation of the 4a-particle condensate state was 
performed with the THSR wave function Q5J . It was conjectured that 4a-particle 
condensation should occur around the 4a threshold. As mentioned in Sec. 12.21 the 
THSR wave function allows only for two limiting configurations, that is a pure 
Slater determinant for B = b and a pure a -particle gas for B^b. Asymptotic con- 
figurations like a+ 12 C(0j~) are absent. Thus, the THSR wave function can well 
describe the dilute a cluster states as well as shell model like ground states, whereas 
other structures such as a+ l2 C(0^) clustering may be only incorporated in an av- 
erage way. In this section, we see whether or not the counterpart of the 0^ state 
obtained by the 4a OCM calculation can also be found with the THSR wave func- 
tion, and then we study how well the 4a condensate state is described with the 
THSR wave function. 
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Table 4 Energies E — E^ a , r.m.s. radii R n ns> monopole matrix elements M(E0), and a decay 
widths r, obtained within the 4a THSR framework, where Ej£v = 4£ K denotes the 4a thresh- 
old energy, with E a the binding energy of the a particle 1221 . 







E-Ef a [MeV] 


#rms [fin] 


M(£0) [fm 2 ] 
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THSR 


-15.05 


2.5 






(0+ 


THSR 


-4.7 


3.1 


9.8 
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THSR 


1.03 


4.2 


2.5 


1.6 


(°4+ 


THSR 


3.04 


6.1 


1.2 


0.14 




The microscopic wave function of ie O with the THSR ansatz is described in 
Eqs. ( TT~8b and dl9) . and the eigenenergies and eigenstates are obtained by solving 
the Hill-Wheeler equation Eq. d20i i. The Hamiltonian is given in Eq. @, where the 
effective nucleon-nucleon interaction called Fl l90l was adopted. The resulting + 
spectrum is shown in Fig. [16] Hereafter, we assign the four0 + states as (O^thsr- 
(O4 )thsr- In Table[H the energies, r.m.s. radii, monopole transition matrix elements 
M(E0) to the ground state, and a-decay widths of the (0|" )thsr-(04~)thsr states are 
displayed. The (O^)thsr state has a large r.m.s. radius of 4.2 fm and the (O^)thsr 
state has an even larger one of 6.1 fm. They are comparable to the values for the 
(O4 )ocm and (0g )ocm states in the 4a OCM calculation, respectively, where the six 
+ states obtained by the 4a OCM calculation in Tableware labeled as (O^ocm- 
(0g )ocm- On the other hand, the M(E0) values of the (O^)thsr and (O4" )thsr states 
well agree with those of the (O4 )ocm and (Of )ocm states, respectively. This sug- 
gests that the (O4 )thsr state corresponds to the (0g )ocm state, and hence to the 
15.1 MeV state. 

More quantitative evidences for the (O4" )thsr state being the counterpart of the 
(0g )ocm state are presented from the analyses of the overlap amplitudes of the a + 
12 C(0 1 h ,0j) channels and the one-body density matrix for the a particle etc. with 
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Fig. 22 (Color online) Radial parts of single-a-particle orbits with L = and = 1 for the 
(O^Jthsr (dotted curve) and (O4 )thsr (solid) states. 



the wave function of the (0^ )thsr state. In Fig.|2T] we show the overlap amplitudes 
of the a + 12 C(0+,0j) channels for the (0j )thsr state. The radial behaviors of the 
a + 12 C(0^,0j) channels are very similar to those of the (O^)ocm case in Fig. IT8l 
In addition, we found that the single-a S orbit occupancy in the (O4 )thsr state is 
as large as 64 %, which is comparable to that in the (Ot )ocm state, and the radial 
behavior of the single-a S orbit in the former state is illustrated in Fig. |22] and is 
similar to that in the latter in Fig. [19] 

In conclusion we found that the (O4 )thsr state corresponds to the (O^)ocm state 
and is most appropriately considered to be the 4a condensate state. This further 
gives us a strong support that the 4a condensate state exists around the 4a breakup 
threshold and is very likely to correspond to the observed 0^ state at 15.1 MeV. 



4.3 Heavier An nuclei: Gross-Pitaevskii Equation 

In principle, one could go on, increasing the number of a-particles, as for 20 Ne, 
24 Mg, etc. and study their structure with use of the THSR wave function or within 
the OCM framework. However, one easily imagines that the complexity of the cal- 
culations quickly becomes prohibitive. In order to get a rough idea what happens for 
more a-particles, drastic approximations have to be performed. One such approxi- 
mation is to consider the a-particles as ideal inert bosons and to treat them in mean 
field approximation. This then leads to the Gross-Pitaevskii Equation (GPE) |9"T1 
which is widely employed in the physics of cold atoms ll83l . One interesting ques- 
tion that can be asked in this connection is: How many a's can maximally exist in a 
self -bound a-gas state? Seeking an answer, it is interesting to investigate it schemat- 
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ically using an effective a-a interaction within an a-gas mean-field calculation of 
the Gross-Pitaevskii type PT1 . 

In the mean field approach, the total wave function of the condensate na-boson 
system is represented as 

<*>(««)= npfr), (49) 

1=1 

where (p is the normalized single-a wave function of the ;-th a boson. Then, the 
equation of motion for the a boson, called the Gross-Pitaevskii equation, is of the 
non-linear Schrodinger-equation type, 

-T- f 1 ~ ~) v2 <P« +^«<P(r) = (50) 
2m \ n J 

C/(r) = (n-l) Jdr'\(p(r')\ 2 V aa (r',r), (51) 

where m stands for the mass of the a particle and U is the mean-field potential 
of a-particles. The center-of-mass kinetic energy correction, 1 — -j, is taken into 
account together with the finite number corrections, n — 1 . In the present study, only 
the 5-wave state is solved self-consistently with the iterative method. The effective 
2a interaction v>2 is taken of Gaussian-type including a repulsive density-dependent 
term, to account for the Pauli repulsion at short distances, which is of similar form 
as the Gogny interaction (known as an effective NN potential) 1921 used in nuclear 
mean-field calculations 

Wr,r')=V exp[-0.7 2 (r-r') 2 ] - 130exp [-0.475 2 (r-r') 2 ] 

+ (47r) 2 gS(r-r')p +V Coo i(r-r'), (52) 

where the units of Vo and r are MeV and fm, respectively, and p denotes the density 
of the na system. The folded Coulomb potential Vcoul is presented as 

Vcoui(r-r') = ^^erf( fl |r-r'|), (53) 
|r-r'| 

The Gaussian-potential part in Eq. ( |52l is based on the Ali-Bodmer poten- 
tial l93l , which is known to reproduce well the elastic a-a scattering phase shift 
up to about 60 MeV for Vo = 500 MeV. The two parameters of the force, Vo and g, 
were adjusted to reproduce the energy (measured from the 3 a threshold) and r.m.s. 
radius of the Hoyle state obtained from the THSR analysis. 

The corresponding a mean-field potential for three a's of 12 C is shown in 
Fig. [23] One sees the OS-state lying slightly above threshold but below the Coulomb 
barrier. As more a-particles are added, the Coulomb repulsion drives the loosely 
bound system of a-particles farther and farther apart. For example, in the case of 
six a's in 24 Mg, the Coulomb barrier is lower and its position is moved outwards. 
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Fig. 23 Alpha-particle mean-field potential for three a's in I2 C and six a's in 24 Mg. The dotted 
line denotes the energy of the single-a 05 orbit (from Ref. HID . 



Thus, eventually the Coulomb barrier should fade away in some limiting nucleus. 
According to our estimate BTI . a maximum of eight to ten a-particles can be held 
together in a condensate. However, there may be ways to lend additional stability 
to such systems. We know that in the case of 8 Be, adding one or two neutrons pro- 
duces extra binding without seriously disturbing the pronounced a-cluster structure. 
Therefore, one has reason to speculate that adding a few neutrons to a many- a state 
may stabilize the condensate. But again, state-of-the-art microscopic investigations 
are necessary before anything definite can be said about how extra neutrons will 
influence an a-particle condensate. 

Concerning excitation of condensate states with many a particles, heavy ion re- 
actions and Coulomb excitation may be appropriate tools. As an ideal case let us 
imagine that 40 Ca has been excited by Coulomb excitation to a state of about 60 
MeV. Coulomb excitation favors + -states and 60 MeV is the threshold for disin- 
tegration into 10 a-particles. Since the Coulomb barrier is absent for ten a's, this 
state may perform a Coulomb explosion of a 10 a particle coherent state. A cartoon 
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of such a scenario is sketched in Fig. [24] With heavy ion reactions, experiments with 
coincident measurements are being analyzed to detect multi a events by Borderie et 
al. (36). W. v. Oertzen et al. seem to have detected an enhancement of multi a decay 
out of a-condensates in compound states of heavier N = Z nuclei, see Ref. |94| . 



4.4 Hoyle-analogue states in non-An nuclei: B and C 

In the previous sections, we discussed the a-gas-like states in An nuclei. On the other 
hand, one can also expect cluster-gas states composed of alpha and triton clusters 
(including valence neutrons etc.) around their cluster disintegrated thresholds in A ^ 
An nuclei, in which all clusters are in their respective OS orbits, similar to the Hoyle 
state with (QS a ) 3 - The states, thus, can be called Hoyle-analogue states in non- 
self-conjugated nuclei. It is an intriguing subject to investigate whether or not the 
Hoyle-analogue states exist in A ^ An nuclei, as for example, B, composed of 2a 
and If clusters as well as 13 C, composed of 3a and In. 

The structure of 3/2" and 1/2+ (1/2 - and 1/2+) states in n B ( 13 C) up to 
around the 2a+t (3a+n) threshold were investigated by the 2a+t OCM |44] 3a+n 
OCM j4~6l ) for 1 *B ( 13 C) combined with the Gaussian expansion method. The model 
space for the 2a+t (3a+«) OCM is large enough to cover the 2a + t (3a + n) gas, 
the 7 Li+a and 8 Be+f ( 9 Be+a and 12 C+«) clusters, as well as the shell-model con- 
figurations. As well known, the a — t and a — n potentials have a strong parity 
dependence J2). In the odd waves they are strongly attractive to produce the bound 
states (for a — t) and resonant states (a — t and a — «), while the even ones are 
weakly attractive and have no ability to produce any resonant states up to E x ~ 15 
MeV |95|. Thus, the gas-like states in n B and 13 C might appear in their even-parity 
states. 

The energy levels of 3/2~ and l/2 + states in n B are shown in Fig. [25] The 
3/2^ state is the ground state with a shell-model-like structure. The calculated nu- 
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Fig. 25 (Color online) Calculated energy levels of 3 /2 and 1/2+ states in 1 1 B with respect to the 
2a + 1 threshold, together with the experimental data. 



clear radius is R = 2.22 fm (R exp = 2.43 ± 0. 1 1 fm). The dominant configuration of 
this state is S\J(3)[f](X,jj,)i = [443] ( 1 , 3) 1 with Q — l harmonic oscillator quanta 
(95 %), having the main angular momentum channel of (L,S)j = (1 , 5)3. On the 

other hand, also the 3 /2^ state has a shell-model-like structure. 

The 3 /2J state appears at E x = 8.2 MeV (E = -2.9 MeV referring to the 2a + 1 
threshold). The radius of 3 /2^ is 3.00 fm. This value is by about 30 % larger than 
that of the ground state of B, and the a — a r.m.s. distance (distance between 
8 Be(2a) and t) is 4.47 fm (3.49 fm). Thus, 3/2^ has a 2a +t cluster structure. A 
characteristic feature of 3 /27 is that the isoscalar monopole transition rate B(E0:IS) 
is as large as 96 ± 16 fm 4 , comparable to that of the Hoyle state (120 ±9 fm 4 ) 8 . The 
present model (92 fm 4 ) reproduces well the data. It is interesting to study whether 
3/2j possesses an a gas nature like the Hoyle state. To this purpose, we study the 
single-cluster orbits and their occupation probabilities in the 3/2^ state by solving 
the eigenvalue equation of the single-cluster density matrices. 

Figure|26]shows the occupation probabilities of the n-th L-wave single- a-particle 
(single-f-particle) orbit in the 3/2^ and 3/27 states. In the 3/2^ state, the occupa- 
tion probabilities of a particles spread out in several orbits, and those of t orbits 
concentrate mainly on two orbits. These results originate from the SU(3) nature of 
the 3/2^ state as mentioned above. On the other hand, in the 3/2^ state, there also 
is no concentration of a occupation probability on a single orbit. This result is in 
contrast with those of the Hoyle state (see Fig.[T2l. Consequently the 3/2^ state can 
not be identified as the analogue of the Hoyle state. The reason why the 3/2^ state 
is not of Hoyle-type is as follows: The 3 /2^ state is bound by 2.9 MeV with respect 
to the 2a + 1 threshold, while the Hoyle state is located by 0.38 MeV above the 3a 
threshold. This extra binding energy of 3/2^ with respect to the 2a + t threshold 
suppresses strongly the development of the gas-like 2a + 1 structure. 
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Fig. 26 (Color online) Occupation probabilities of the a (?) orbits for the (a) 3/2j and (b) 3/2 3 
states in "B. 



As for the 1 /2+ states, the 1 /2~\ state appears as a bound state at £f p = 6.79 
MeV around the 7 Li+a threshold. This low excitation energy indicates that a -type 
correlations should play an important role in this state. In fact, we found that the 
1 /2j" state with the radius of 3.14 fm has a 7 Li(g.s)+a structure with P-wave rela- 
tive motion, although the Li(a + 1) part is rather distorted in comparison with the 
ground state of 7 Li. Since the 3/2^ state has the largest S" factor for the 7 Li(g.s)+oc 
channel with 5-wave relative motion, the 1/2^ and 3/27 states of B can be inter- 
preted as parity-doublet partners. 

In addition to the 1/2+ state, the 1/2+ state appears as a resonant state at 
E x = 11.95 MeV (r = 190 keV) around the 2a + t threshold using the complex- 
scaling method 196TI991 . The large radius (Rn = 5.98 fm) indicates that the state has 
a dilute cluster structure. The analysis of the single-cluster properties showed that 
this state has as main configuration (0S a ) 2 {0S t ) orbital occupation with about 65% 
probability (see Fig. |26V Thus, the 1/2+ state can be called the Hoyle-analogue. 
Recently, the 1/2+ (3/2+) state at E x = 12.56 MeV with T = 210 ±20 keV (lo- 
cated at 1.4 MeV above the 2a + t threshold) was observed in the oc+ 7 Li decay 
channel [ 100-102 1 . The energy and width of the 12.56-MeV state are in good cor- 
respondence to the present study. The Hoyle-analogue state in B, thus, could be 
assigned as the 12.56-MeV state. It should be reminded that the 1/2+ state is lo- 
cated by 0.75 MeV above the a + a + 1 threshold, while 1 /2 + is bound by 4.2 MeV 
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Fig. 27 (Color online) Calculated energy levels of 1 /2 and 1/2+ states in 13 C with respect to the 
3a + f threshold, together with the experimental data. 



with respect to the three cluster threshold. The latter binding energy leads to a sup- 
pression of the development of the gas-like a + a + t structure in 1 j2\, whereas it 
is generated with a large nuclear radius in the 1 j2\ state because of its appearance 
above the three-body threshold. 

The calculated energy spectrum of 1 /2~ states in 13 C is shown in Fig. [27] The 
four 1 /2~ energy levels are in good correspondence with the experimental data. The 
ground state (l/2p is described as having a shell-model configuration. The calcu- 
lated nuclear radius (2.39 fm) agrees with the data (2.44 fm). The three isoscalar 
monopole transition strengths, M(l/2J - 1/2^ ) = 4.2, M{\j2\ - l/2^~) = 5.6, 
and M(l/2^ — 1/2^) = 8.2, are also consistent with experiment 0451 . 6.1 ±0.5, 
4.2 ± 0.4, and 4.9 ± 0.4, respectively, in units of fm 2 . The nuclear radii for the three 
excited 1 /2~ states are 3.36, 2.96, and 3.19 fm, respectively. 

From the analysis of the radial behavior of the overlap amplitudes referring to 
the 12 C+« and 9 Be+oc channels, the 1 /2^ and 1 /2^ states are characterized as hav- 
ing large components of I2 C(g.s,2 + )+n and 12 C(3~)+n with 9 Be(g.s)+ot, respec- 
tively. On the other hand, the 1 /2^ state contains a somewhat large component of 
the 12 C(Hoyle)+« channel together with 12 C(2 + )+« and 9 Be(g.s)+a. However, this 
state does not have as large an a condensate component as the Hoyle state in 12 C. 
This is due to strong attraction of the odd-wave a — n potentials which induces the 
coupling of 12 C(2 + )+« and 9 Be(g.s)+a structures with the 12 C(Hoyle)+« configu- 
ration, and disturbs significantly the structure of the 3 a condensate in 13 C. In the 
mirror nucleus 13 N, the 3a+p OCM analysis gives qualitatively similar results to 
those of the present 13 C case 11031 . 

As for the 1/2+ states of 13 C, the energy spectrum of the first three 1 /2 + states 
correspond well with the data (see Fig. l27t . We found that the l/2j~ state has a 
main configuration [ 12 C(g.s) (g> 2si^]- Reflecting the fact that the neutron binding 
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energy of the 1/2|~ state with respect to the 12 C(g.s.)+n threshold is as small as 
1.9 MeV, this state has a neutron-halo-like structure. In fact the calculated nuclear 
radius of this state (2.68 fm) is larger than that of the ground state (2.39 fm), and 
this enhancement of the radius comes from the neutron-halo- like structure. 

On the other hand, the 1 /2t state has a dominant configuration of the extra neu- 
tron coupled with the Hoyle state, with non-negligible mixing of 9 Be(g.s,l /2j~)+a 
channels. The nuclear radius is about 4.0 fm, which is smaller than that of the Hoyle 
state in the 3a OCM calculation lETTl . We found that the size of the 3a part in this 
state is reduced by about 15 % in comparison with that of the Hoyle state. The oc- 
cupation probability of a particle in OS orbit in this state is less than 30 %, which is 
much smaller than that for the Hoyle state. 

The 1 state around the 3<x + n threshold has the nuclear radius of 5.40 fm with 
a dilute a condensate feature, in which 3 a particles occupy an identical OS orbit 
with 55 % probability. This state has a rather large overlap with the 9 Be(l /2 1 h )+a 
channel as well as with the 12 C(Hoyle)+« one. It is noted that the 9 Be(l/2 1 h ) state 
is known to have a neutron-halo-like structure (or 2a + n gas-like structure). Thus, 
these results suggest that the 1/2^ state is a candidate for the Hoyle-analogue state. 

With this we terminate our consideration of cluster and condensate aspects in 
finite nuclei. An important connection with the finite systems is given by clustering, 
for instance a particle clustering and condensation in infinite matter. In the next 
section we turn to these issues. 



5 Clusters in nuclear matter and a -particle condensation 
5.1 Nuclear clusters in the medium 

Of course, it is also interesting and important to study how a -clusters behave and 
actually condensate in infinite symmetric and asymmetric nuclear matters. This not 
only in regard to better understand what finally happens in a finite nucleus but in 
collapsing and compact stars one may speculate about the existence of a macro- 
scopic a-particle condensate. So let us first consider the modification an a particle 
undergoes when it is embedded in a nuclear medium. 

Medium modifications of single-particle states as well as of few-nucleon states 
become of importance with increasing density of nuclear matter. The self-energy 
of an A-particle cluster can in principle be deduced from contributions describing 
the single-particle self-energies as well as medium modifications of the interac- 
tion and the vertices. A guiding principle in incorporating medium effects is the 
construction of consistent ("conserving") approximations, which treat medium cor- 
rections in the self-energy and in the interaction vertex at the same level of ac- 
curacy. This can be achieved in a systematic way using the Green functions for- 
malism 111041 . At the mean-field level, we have only the Hartree-Fock self-energy 
F HF (1) = L2^(12, 12)/(2) together with the Pauli blocking factors, which modify 
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the two-nucleon interaction from V(12,1'2') to V(12,l'2')[l -/(l) ~/(2)], with 
/(l) = [1 +exp(£ HF (l) -n)/T]- x and V{\2, 12) = V(12, 12) - V(12,21). In the 
case of the two-nucleon system (A = 2), the effective wave equation which includes 
those corrections is presented in the following form \ 1051 11 06l . 

[e HF (l) + e HF (2)-£ 2 ,p] v/ 2 .p(12) 

+ i£[l-/(l)-/(2)]V(12,l'2')w(l'2')=0. (54) 

Z 1'2' 

This effective wave equation describes bound states as well as scattering states. The 
onset of pair condensation is achieved when the binding energy £'2,/>=o coincides 
with 2/1, where P denotes the total momentum of the two-nucleon system. It is 
noted that the Gor'kov equation in BCS theory of superfluidity is a special case of 
Eq. dH. 

Similar equations have been derived from the Green function approach for the 
case A = 3 and A = 4, describing triton/helion ( 3 He) nuclei as well as a -particles in 
nuclear matter. The effective wave equation contains in mean field approximation 
the Hartree-Fock self-energy shift of the single-particle energies as well as the Pauli 
blocking of the interaction. We give the effective wave equation for A = 4, 

[e HF (l) + e HF (2) + e HF (3) + e HF (4) -E Ar ] yM 1234 ) 

+il E [i-/(o-/(j)]n'j,'7)nvft,p(iw)=o. (55) 

Z /</l'2'3'4' k£i,j 

A similar equation is obtained for A = 3, which is an equation for a fermionic cluster. 

The effective wave equation has been solved using separable potentials for A = 2 
by integration. For A = 3 , 4 we can use a Faddeev approach ||49| . The shifts of bind- 
ing energy can also be calculated approximately via perturbation theory. In Fig. [28] 
we show the shift of the binding energy of the light clusters (d,t/h and a) in sym- 
metric nuclear matter as a function of density for temperature T = 10 MeV |49l . 

It is found that the cluster binding energy decreases with increasing density. Fi- 
nally, at the Mott density p^p tl (T) the bound state is dissolved. The clusters are not 
present at higher densities, merging into the nucleonic medium. For a given cluster 
type characterized by A,«, we can also introduce the Mott momentum P^ ott (p,T) 
in terms of the ambient temperature T and nucleon density p, such that the bound 
states exist only for P > P^ lott (p , T). We do not present an example here, but it is in- 
tuitively clear that a cluster with high c.o.m. momentum with respect to the medium 
is less affected by the Pauli principle than a cluster at rest, because the overlap of the 
bound state wave function in momentum space and the Fermi distribution function 
becomes smaller. 
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Fig. 28 Shift of binding energy of the light clusters (d - dash dotted, t/h - dotted, and a - dashed: 
perturbation theory, full line: non-perturbative Faddeev-Yakubovski equation) in symmetric nu- 
clear matter as a function of density for given temperature T = 10 MeV 1491 . 



5.2 Four-particle condensates and quartetting in nuclear matter 

In general, it is necessary to take into account all bosonic clusters to gain a com- 
plete picture of the onset of superfluidity. As is well known, the deuteron is weakly 
bound as compared to other nuclei. Higher A-clusters can arise that are more stable. 
In this section, we will consider the formation of a -particles, which are of special 
importance because of their large binding energy per nucleon (~ 7 MeV). We will 
not include tritons or helions, which are fermions and not so tightly bound. More- 
over, we will not consider nuclei in the iron region, which have even larger binding 
energy per nucleon than the a -particle and thus constitute, in principle, the domi- 
nant component at low temperatures and densities. However, the latter are complex 
structures of many particles and are strongly affected by the medium as the density 
increases for given temperature, so that they are assumed not to be of relevance in 
the density region considered here. 

The in-medium wave equation for the four-nucleon problem has been solved 
using the Faddeev-Yakubovski technique, with the inclusion of Pauli blocking, see 
also below. The binding energy of an a-like cluster with zero c.o.m. momentum van- 
ishes at around po/10, where po ~ 0.16 nucleons/fm 3 denotes the saturation density 
of isospin-symmetric nuclear matter, see Fig. |28] Thus, the four-body bound states 
make no significant contribution to the composition of the system above this density. 
Given the medium-modified bound-state energy £4 p, the bound-state contribution 
to the EOS is 



We will not include the contribution of the excited states nor that of scattering states. 
Because of the large specific binding energy of the a particle, low-density nuclear 



1 -1 



(56) 
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matter is predominantly composed of a particles. This observation underlies the 
concept of a matter and its relevance to diverse nuclear phenomena H107H1 1 ll . 

As exemplified by Eq. d55l l. the effect of the medium on the properties of an 
a particle in mean-field approximation (i.e., for an uncorrelated medium) is pro- 
duced by the Hartree-Fock self-energy shift and Pauli blocking. The shift of the 
a-like bound state has been calculated using perturbation theory [105 106] as well 
as by solution of the Faddeev-Yakubovski equation l49l . It is found that the bound 
states of clusters d, t, and h with A < 4 are already dissolved at a Mott density 
pMott p /io, see Fig. [28] Since Bose condensation only is of relevance for d and 
a, and the fraction of d, t and h becomes low compared with that of a with in- 
creasing density, we can neglect the contribution of them to an equation of state. 
Consequently, if we further neglect the contribution of the four-particle scattering 
phase shifts in the different channels, we can now construct an equation of state 
p(T,n) = p bee {T,ii) + p boundM (T,ii)+p hoand - a (T,ii) such that a-particles deter- 
mine the behavior of symmetric nuclear matter at densities below p%} ott and tem- 
peratures below the binding energy per nucleon of the a-particle. The formation of 
deuteron clusters alone gives an incorrect description because the deuteron binding 
energy is small, and, thus, the abundance of cZ-clusters is small compared with that 
of a-clusters. In the low density region of the phase diagram, a-matter emerges as 
an adequate model for describing the nuclear-matter equation of state. 

With increasing density, the medium modifications - especially Pauli blocking - 
will lead to a deviation of the critical temperature T c (p ) from that of an ideal Bose 
gas of a-particles (the analogous situation holds for deuteron clusters, i.e., in the 
isospin-singlet channel) l49l . 

Symmetric nuclear matter is characterized by the equality of the proton and neu- 
tron chemical potentials, i.e., fi p = fi„ = fi . Then an extended Thouless condition 
based on the relation for the four-body T-matrix (in principle equivalent to Eq. (T55t 
at eigenvalue 4-ji) 



serves to determine the onset of Bose condensation of a-like clusters, noting that 
the existence of a solution of this relation signals a divergence of the four-particle 
correlation function. An approximate solution has been obtained by a variational 
approach, in which the wave function is taken as Gaussian incorporating the correct 
solution for the two-particle problem 048] ■ 

On the other hand, Eq. d57l l. respectively Eq. (l55T l at eigenvalue 4/x, has also 
been solved numerically exactly by the Faddeev- Yakubovsky method employing the 
Malfliet-Tjon force Ml 121 . The results for the critical temperature of a-condensation 
is presented in Fig.|29]as a function of the chemical potential ji (see also Ref. l48l ). 
The exact solution could only be obtained for negative /i, i.e. when there exists a 



T 4 (1234,l"2"3"4",4jLi) = \ 




V(12,l'2')[l-/(l)-/(2)] 
4 J u-£ , i-£ , 2 -£ , 3-£ , 4 



5(3,305(4,4') 
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Fig. 29 Critical temperature of alpha and deuteron condensations as functions of (a) chemical 
potential and (b) density of free nucleons |50|. Crosses (x) correspond to the solution of Eq. j55t 
with the Malfliet-Tjon interaction (MT I-III) using the Faddeev-Yakubovski method. 



bound cluster. It is, therefore, important to try yet another approximate solution of 
the in-medium four-body equation. Since the a-particle is strongly bound, we make 
a momentum projected mean field ansatz for the quartet wave function 1113141151 

4 

^1234 = (2^) 3 5^(k! +k 2 + k 3 +k 4 )n'P(k,)Z Sr , (58) 

1=1 

where % ST is the spin-isospin function which we suppose to be the one of a scalar 
(S = T = 0). We will not further mention it from now on. We work in momentum 
space and <p(k) is the as-yet unknown single particle OS wave function. In position 
space, this leads to the usual formula |fl9l fi234 — > J ^ 3 ^Ilf=i^( r ; — R) where 
<p(r,) is the Fourier transform of <p(k,). If we take for <p(k,-) a Gaussian shape, 
this gives: fi234 — > ^Vp[~ c Hi<i<k<A{ r i ~ r A-) 2 ] which is the translationally invariant 
ansatz often used to describe a-clusters in nuclei. For instance, it is also employed 
in the a-particle condensate wave function of Tohsaki, Horiuchi, Schuck, Ropke 
(THSR)inRef. 0. 

Inserting the ansatz (l58l into (T55t and integrating over superfluous variables, or 
minimizing the energy, we arrive at a Hartree-Fock type of equation for the single 
particle OS wave function <p (k) = <p (|k|) which can be solved. However, for a general 
two body force Vj^k, t' k' > tne equation to be solved is still rather complicated. We, 
therefore, proceed to the last simplification and replace the two body force by a 
unique separable one, that is 
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V k|k ,^ k ^A^ 2 ^^' 2 /^ (2 ^)3 5 (3) (K _K'), (59) 



where k = (kj - k 2 )/2, k' = (k'j - k 2 )/2, K = kj +k 2 , and K' = k[ +k' 2 . This 
means that we take a spin-isospin averaged two body interaction and disregard that 
in principle the force may be somewhat different in the S, T = 0, 1 or 1 , channels. 
It is important to remark that for a mean field solution the interaction only can be 
an effective one, very different from a bare nucleon-nucleon force. This is contrary 
to the usual gap equation for pairs, to be considered below, where, at least in the 
nuclear context, a bare force can be used as a reasonable first approximation. 

We are now ready to study the solution of Eq. (I5"5l l for the critical temperature 
T c a , defined by the point where the eigenvalue equals 4fi. For later comparison, 
the deuteron (pair) wave function at the critical temperature is also deduced from 
Eqs. (|55]) and (|59]l to be 



(60) 



where (k) is the relative wave function of two particles given by ffo — > <j> ( | 2 I ) 
5' 3 '(kj + k 2 ), and e = k 2 / (2m). We also neglected the momentum dependence of 
the Hartree-Fock mean field shift in Eq. (|60V It, therefore, can be incorporated into 
the chemical potential jj. . With Eq. d60l i. the critical temperature of pair condensation 
is obtained from the following equation: 

J (27T) 3 k 2 /m-2ix 

In order to determine the critical temperature for a-particle condensation, we 
have to adjust the temperature so that the eigenvalue of (IB3T l and d57l i equals 4-jj.. 
The result is shown in Fig. |29f a). In order to get an idea how this converts into a 
density dependence, we use for the moment the free gas relation between the density 
of uncorrected nucleons and the chemical potential 

We are well aware of the fact that this is a relatively gross simplification, for in- 
stance at the lowest densities, and we intend to generalize our theory in the future 
so that correlations are included into the density. This may be done along the work 
of Nozieres and Schmitt-Rink 111 161 . The two open constants A and ko in Eq. d59l 
are determined so that binding energy (—28.3 MeV) and radius (1.71 fm) of the free 
(fi = 0) a-particle come out right. The adjusted parameter values are: X = —992 
MeV fm 3 , and ko = 1.43 fm -1 . The results of the calculation are shown in Fig. [29] 
In Fig. [29] the maximum of critical temperature T c a max is at fi = 5.5 MeV, and the 
a-condensation can exist up to /i max = 11 MeV. It is very remarkable that the results 
obtained with d58l for T t a very well agree with the exact solution of (I5"5l l and (IBTt 
using the Malfiiet-Tjon interaction (MT I-III) 111 121 with the Faddeev-Yakubovski 
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Fig. 30 Single particle wave functions (al~a4) in momentum space (p(k) and (bl~b4) in position 
space r(p(r) at chemical potential (/I), critical temperature (T c ), and density (n), which are obtained 
by solving Eq. {55} with the mean field ansatz jSjjQ ESI: for (al) [(bl)] jX = -7.08 MeV, T c = 
MeV, n = fm" 3 , for (a2) [(b2)] /I = -2.22 MeV, T c = 6.61 MeV, n = 9.41 x 1(T 3 fm" 3 , for (a3) 
[(b3)] 11 = 6.11 MeV, F t . = 8.45 MeV, n = 3.07 x 10~ 2 fm" 3 , and for (a4) [(b4)] li = 10.6 MeV, 
T c = 5.54 MeV, n = 3.34 x 10~ 2 fm~ 3 . Figs, (al) and (bl) correspond to the wave functions for 
free a-particle. The vertical lines in Figs. (a3) and (a4) are at the Fermi wave length kf = ^/2m\i. 



method also shown by crosses in Fig. [29] (the numerical solution only could be 
obtained for negative values of fi). This indicates that T t a is essentially determined 
by the Pauli blocking factors. 

In Fig. [29] we also show the critical temperature for deuteron condensation de- 
rived from Eq. (|6TV In this case, the bare force is adjusted with A = —1305 MeV 
fm 3 and Icq = 1.46 fm -1 to get experimental energy (—2.2 MeV) and radius (1.95 
fm) of the deuteron. It is seen that at higher densities deuteron condensation wins 
over the one of a-particle. The latter breaks down rather abruptly at a critical posi- 
tive value of the chemical potential. Roughly speaking, this corresponds to the point 
where the a-particles start to overlap. This behavior stems from the fact that Fermi- 
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Dirac distributions in the four body case, see Eq. (l55l , can never become step-like, 
as in the two body case, even not at zero temperature, since the pairs in an a-particle 
are always in motion. As a consequence, a-condensation generally only exists as a 
BEC phase and the weak coupling regime is absent, see also discussion in Sec. 15.41 

Figure [30] shows the normalized self-consistent solution of the wave function in 
momentum space derived from Eq. (IB3T l with the mean field ansatz ( f58l ) and the wave 
function in position space defined by its Fourier transform <p{r). Figures [30t al) and 
(bl) represent the wave functions of the free a-particle. The wave function resem- 
bles a Gaussian and this shape is approximately maintained as long as ji is negative, 
see Fig. [30f a2). On the contrary, the wave function of Fig. I307a3). where the chem- 
ical potential is positive, has a dip around k = which is due to the Pauli blocking 
effect. For the even larger positive chemical potential of Fig. [3Ql a4) the wave func- 
tion develops a node. The maximum of the wave function shifts to higher momenta 
and follows the increase of the Fermi momentum kf, as indicated on Fig. [30] 

On the other hand, the wave functions in position space in Figs. I30l b2). (b3) 
and (b4) develop an oscillatory behavior, as the chemical potential increases. This 
is reminiscent to what happens in BCS theory for the pair wave function in position 
space 11171 . 

An important consequence of this study is that at the lowest temperatures, Bose- 
Einstein condensation occurs for a particles rather than for deuterons. As the den- 
sity increases within the low-temperature regime, the chemical potential p first 
reaches —7 MeV, where the a's Bose-condense. By contrast, Bose condensation 
of deuterons would not occur until ji rises to —1.1 MeV. 

The " quartetting " transition temperature sharply drops as the rising density ap- 
proaches the critical Mott value at which the four-body bound states disappear. At 
that point, pair formation in the isospin-singlet deuteron-like channel comes into 
play, and a deuteron condensate will exist below the critical temperature for BCS 
pairing up to densities above the nuclear-matter saturation density po, as described 
in the previous Section. Of course, also isovector n-n and p-p pairing develops. The 
critical density at which the a condensate disappears is estimated to be po/3. There- 
fore, a-particle condensation primarily only exists in the Bose-Einstein-Condensed 
(BEC) phase and there does not seem to exist a phase where the quartets acquire a 
large extension as Cooper pairs do in the weak coupling regime. However, the vari- 
ational approaches of Ref. l48l and of Eq. (l58l l on which this conclusion is based 
represent only a first attempt at the description of the transition from quartetting 
to pairing. The detailed nature of this fascinating transition remains to be clarified. 
Many different questions arise in relation to the possible physical occurrence and 
experimental manifestations of quartetting: Can we observe the hypothetical "a 
condensate" in nature? What about thermodynamic stability? What happens with 
quartetting in asymmetric nuclear matter? Are more complex quantum condensates 
possible? What is their relevance for finite nuclei? As discussed, the special type 
of microscopic quantum correlations associated with quartetting may be important 
in nuclei, its role in these finite inhomogeneous systems being similar to that of 
pairing. 
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Fig. 31 Critical temperature as a function of the total chemical potential /i t otal = M;> + Mn (top) 
and the total free density )! to tal (bottom) [51 1. Thick (thin) lines are for a -particle (deuteron). Solid, 
dashed, and dotted lines are respectively for 5 = 0.0, 5 = 0.5, and 5 = 0.9, where the density ratio 
5 is as in Eq. )64l i. 



On the other hand, if at all, a-condensation in compact star occurs at strongly 
asymmetric matter. It is, therefore, important to generalize the above study for 
symmetric nuclear matter to the asymmetric case. This can be done straight for- 
wardly again using our momentum projected mean field ansatz (l58l l generalized to 
the asymmetric case. This implies to introduce two chemical potentials, one for neu- 
trons and for protons. We also have to distinguish two single particle wave functions 
in our product ansatz which now reads 

Vi234 -> (Pp(ki)(p p {k 2 )(p n (k 3 )q) n (k4)xo 

x (27r) 3 5(k 1 +k 2 + k 3 + k 4 ) (63) 

where <p T (k,) = <p T (|k,|) is the s-wave single particle wave functions for protons (t = 
p) and neutrons (t = n), respectively. Xo is the spin-isospin singlet wave function. 
This now leads to two coupled equations of the Hartree-Fock type for <p„ and p„. 
For the force we use the same as in the symmetric case. 

Fig- EH a ) shows the critical temperature of a condensation as a function of the 
total chemical potential /i to tai = + /V We see that T c decreases as the asymmetry, 



52 



T. Yamada, Y. Funaki, H. Horiuchi, G. Ropke, P. Schuck, and A. Tohsaki 



given by the parameter 



5 



n 



n 



p 



n„ + n 



(64) 



p 



increases. This is in analogy with the deuteron case (also shown) which already had 
been treated in Refs. 11 1811 191 . On the other hand, in Fig. l3TT b). it is also interesting 
to show T c as a function of the free density which is 



where the factor two in front of the integral comes from the spin degeneracy, and 
fp,n(k) = [1 + exp(7i 2 £ 2 /2ra — jttp. n )] _1 . It should be emphasized, however, that in the 
above relation between density and chemical potential, the free gas relation is used 
and correlations in the density have been neglected. In this sense the dependence of 
T c on density only is indicative, more valid at the higher density side. The very low 
density part where the correlations play a more important role shall be treated in a 
future publication. It should, however, be stressed that the dependence of T c on the 
chemical potential as in Fig. l3TT a). stays unaltered. 

The fact that for more asymmetric matter the transition temperature decreases, 
is natural, since as the Fermi levels become more and more unequal, the proton- 
neutron correlations will be suppressed. For small 5's, i.e., close to the symmetric 
case, a condensation (quartetting) breaks down at smaller density (smaller chem- 
ical potential) than deuteron condensation (pairing). This effect has already been 
discussed in our previous work for symmetric nuclear matter [48 50|. For large 5's, 
i.e. strong asymmetries, the behavior is opposite, i.e., deuteron condensation breaks 
down at smaller densities than a condensation, because the small binding energy of 
the deuteron can not compensate the difference of the chemical potentials. 

More precisely, for small 5's, the deuteron with zero center of mass momentum 
is only weakly influenced by the density or the total chemical potential as can seen in 
Fig. [3T| However, as 8 increases, the different chemical potentials for protons and 
neutrons very much hinders the formation of proton-neutron Cooper pairs in the 
isoscalar channel for rather obvious reasons. The point to make here is that because 
of the much stronger binding per particle of the a-particle, the latter is much less 
influenced by the increasing difference of the chemical potentials. For the strong 
asymmetry 8 = 0.9 in Fig. [3T1 then finally a-particle condensation can exist up to 
itotai = 0.02 fm~ 3 (/i to tai = 9.3 MeV), while the deuteron condensation exists only 
up to n total = 0.005 far 3 (itotai = 6.0 MeV). 

Overall, the behavior of T c is more or less as expected. We should, however, 
remark that the critical temperature for a-particle condensation stays quite high, 
even for the strongest asymmetry considered here, namely 8 = 0.9. This may be 
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Fig. 32 Momentum-space single particle wave functions for proton <p p (solid line), and for neutron 
<p„ (dashed line) for the critical temperature as a function of k for 5 = 0.0, 0.5 0.9, and the real- 
space wave functions for proton <p p (solid line), for neutron <p„ (dashed line) as a function of r for 
5 = 0.9 derived from the Fourier transform of (p pj , (k) with <p pn (r) = / d 3 ke' k ' T (p p „ (k) / {2k) 3 BTI . 
The top, middle and bottom figures are for ^i t otal = lh> + Vhi ~ — 1 1 MeV, ~ 0.0 MeV, and ~ 
9.0MeV, respectively. The wave functions are normalized by / cPk(p^ n (k) / '(2tt) 3 = 1. The details 
of data for respective figures are following: (al) 8 = 0.0, /X to tal = — 1 1 • 1 MeV, [l p = —5.53 MeV, 
\ln = -5.53 MeV, T c = 4.52 MeV. (a2) 8 = 0.5, /i tota i = -11.5 MeV, jx p = -8.18 MeV, 
H„ = -3.35 MeV, T c = 4.07 MeV. (a3), (a4) 8 = 0.9, /i tota i = -11.0 MeV, \l p = -10.8 MeV, 
ll n = -0.163 MeV, T L . = 3.35 MeV. (bl) 5 = 0.0, ^total = 0.028 MeV, \l p = -0.014 MeV, \i n = 
-0.014 MeV, T c = 7.46 MeV. (b2) 8 = 0.5, ^otal = 0.11 MeV, ll p = -4.65 MeV, fl„ = 4.76 MeV, 
T c = 6.74 MeV. (b3), (b4) 5 = 0.9, ^, ota i = -0.02 MeV, jl p = -8.18 MeV, ^„ = 8.16 MeV, 
T c = 4.29 MeV. (cl) 8 = 0.0, ^otal = 8.80 MeV, p. p = 4.40 MeV, ^„ = 4.40 MeV, T c = 8.44 MeV. 
(c2) 8 = 0.5, /i,otal = 8.93 MeV, ^ = -1.12 MeV, \i n = 10.0 MeV, T, = 7.16 MeV. (c3), (c4) 
5 = 0.9, /i tota i = 8.94 MeV, \i p = -4.21 MeV, \l n = 13.2 MeV, T c = 3.72 MeV. 



of importance for the possibility of a-particle condensation in neutron stars and 
supernovae explosions H20l ll2ll . 

We also show the single particle wave functions of protons and neutrons, en- 
tering the quartet wave function (163) . for various ratios of Fermi surface imbal- 
ance and chemical potentials in Fig. [32] In most cases of Fig. [32J the momentum- 
space wave functions with negative chemical potentials are monotonically decreas- 
ing whereas the ones with positive chemical potentials have a dip at k = 0. However, 
the momentum-space wave functions also develop a dip at k = even at a negative 
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chemical potential as the asymmetry takes on stronger values. This can be seen in 
Fig. |32f a3) and (c2). Furthermore, the neutron wave function in £-space with large 
positive chemical potential develops a node. This behavior is similar to the wave 
functions in Ref. l50l . As shown in Fig.[32l the dissymmetry of proton and neutron 
wave functions increases as 8 increases. As a consequence, the critical temperature 
decreases, and the a condensation breaks down at a more dilute density, see Fig.l3Tl 
We also present in Fig. [321" a4). (b4) and (c4) the proton and neutron wave functions 
in real space. In spite of the sometimes strong dissymmetry in momentum space, 
the proton and neutron wave functions are relatively more similar to one another in 
r-space. The neutron wave function develops a node as the total chemical potential 
Mtotai = + M» increases, but the negative values of the wave function remain rather 
moderate. 

In conclusion the a-particle (quartet) condensation was investigated in homoge- 
neous symmetric nuclear matter as well as in asymmetric nuclear matter. We found 
that the critical density at which the a-particle condensate appears is estimated to 
be around Po/3 in the symmetric nuclear matter, and the a-particle condensation 
can occur only at low density. This result is consistent with the fact that the Hoyle 
state (Oj) of I2 C also has a very low density p ~ po/3. On the other hand, in the 
asymmetric nuclear matter, the critical temperature T c for the a-particle conden- 
sation was found to decrease with increasing asymmetry. However, T c stays rela- 
tively high for very strong asymmetries, a fact of importance in the astrophysical 
context. The asymmetry affects deuteron pairing more strongly than a-particle con- 
densation. Therefore, at high asymmetries, if at all, a-particle condensate seems to 
dominate over pairing at all possible densities. 



5.3 Reduction of the a-condensate with increasing density 

The properties of a matter can be used to frame the discussion of the structure of 
n a nuclei. As described in the preceding section, computational studies of these 
nuclei based on THSR cluster states have demonstrated that an a condensate is es- 
tablished at low nucleon density. More specifically, states lying near the threshold 
for decomposition into a particles, notably the ground state of 8 Be, 12 C in the 0^ 
Hoyle state, and corresponding states in ie O and other na nuclei are dilute, being of 
low mean density and unusually extended for their mass numbers. We have shown 
quantitatively within a variational approach that a-like clusters are well formed, 
with the pair correlation function of a-like clusters predicting relatively large mean 
distances. For example, in determining the sizes of the 12 C nucleus in its 0^ (ground) 
state and in its 0^ excited state, we obtained the r.m.s. radii of 2.44 fm and 3.83 fm, 
respectively. The corresponding mean nucleon densities estimated from 36/4?rr 3 ms 
are close to the nuclear-matter saturation density po = 0.16 nucleon/fm 3 in the for- 
mer state and 0.03 nucleon/fm 3 in the latter. The expected low densities of putative 
alpha-condensate states are confirmed by experimental measurements of form fac- 
tors 022). 
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All of our considerations indicate that quartetting is possible in the low-density 
regime of nucleonic matter, and that a condensates can survive until densities of 
about 0.03 nucleons/fm 3 are reached. Here, we are in the region where the concept 
of a matter can reasonably be applied 0123111241 . It is then clearly of interest to 
use this model to gain further insights into the formation of the condensate, and 
especially the reduction or suppression of the condensate due to repulsive interac- 
tions [88]. We will show explicitly that in the model of a matter, as in our studies of 
finite nuclei, condensate formation is diminished with increasing density. Already 
within an a-matter model based on a simple a — a interaction, we can demon- 
strate that the condensate fraction - the fraction of particles in the condensate - is 
significantly reduced from unity at a density of 0.03 nucleon/fm 3 and essentially 
disappears approaching nuclear matter-saturation density. 

The quantum condensate formed by a homogeneous interacting boson system at 
zero temperature has been investigated in the classic 1956 paper of Penrose and On- 
sager 11251 who characterize the phenomenon in terms of off-diagonal long-range 
order of the density matrix. Here we recall some of their results that are most rele- 
vant to our problem. Asymptotically, i.e., for |r — r'| ~°o, the nondiagonal density 
matrix in coordinate representation can be decomposed as 

p (r, r') ~ y/ * (r)^o(r') + y(r - r') . (68) 

In the limit, the second contribution on the right vanishes, and the first approaches 
the condensate fraction, formally defined by 

Penrose and Onsager showed that in the case of a hard-core repulsion, the con- 
densate fraction is determined by a filling factor describing the ratio of the volume 
occupied by the hard spheres. They applied the theory to liquid 4 He, and found that 
for a hard-sphere model of the atom-atom interaction yielding a filling factor of 
about 28%, the condensate fraction at zero temperature is reduced from unity (its 
value for the noninteracting system) to around 8%. (Remarkably, but to some extent 
fortuitously, this estimate is in rather good agreement with current experimental and 
theoretical values for the condensate fraction in liquid 4 He.) 

To make a similar estimate of the condensate fraction for a matter, we follow 
Ref. 11211 and assume an "excluded volume" for a particles of 20 fm 3 . At a nucle- 
onic density of po/3, this corresponds to a filling factor of about 28%, the same as 
for liquid 4 He. Thus, a substantial reduction of the condensate fraction from unity 
(for a noninteracting a-particle gas at zero temperature) is also expected in low- 
density a matter. 

Turning to a more systematic treatment, we proceed in much the same way 
as Clark and coworkers 01231 . referring especially to the most recent study with 
M. T. Johnson. Adopting the a — a interaction potential 

Vaa(r) = 475 e - (0 - 7r/fm)2 MeV - 130 e -(°- 475r/fm > 2 MeV (70) 
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introduced by Ali and Bodmer 11931 , we calculate the reduction of the condensate 
fraction as function of density within what is now a rather standard variational ap- 
proach. Alpha matter is described as an extended, uniform Bose system of interact- 
ing a particles, disregarding any change of the internal structure of the a clusters 
with increasing density. In particular, the dissolution of bound states associated with 
Pauli blocking (Mott effect) is not taken into account in the present description. 

The simplest form of trial wave function incorporating the strong spatial correla- 
tions implied by the interaction potential d70t is the familiar Jastrow choice, 

^(r 1 ,...,r A )=n/(|r ! --r 7 -|). (71) 

The normalization condition 

4np a r\f 2 {r)-\]r 2 dr = -\, (72) 
Jo 

in which p a is the number density of a-particles, is imposed as a constraint on 
the variational wave function, in order to promote the convergence of the cluster 
expansion used to calculate the energy expectation value [ 1 26] . In the low-density 
limit, the energy functional [binding energy per a cluster as a functional of the 
correlation factor f(r)] is given by 

E[f]=27Z Pa £^ (^My + f(r)V a (r)Ydr, (73) 

where m a is the a-particle mass, while the condensate fraction is given by 

Po = exp^-4n Pa J~[f{r)-l} 2 r 2 drj . (74) 

The variational two-body correlation factor / was taken as one of the forms em- 
ployed by Clark and coworkers 01231 , namely 

f(r) = (1 -e- c "-)(l+be- ar + ce- 2ar ). (75) 

At given density p, the expression for the energy expectation value is minimized 
with respect to the parameters a, b, and c, subject to the constraint d72l . It is impor- 
tant to note that these approximations, based on truncated cluster expansions, are 
reliable only at densities low enough so that the length scale associated with decay 
of f 2 ( 1 ) — 1 is sufficiently small compared to the average particle separation, which 
is inversely proportional to the cubic root of the density [123, 124 1261 11271 . 

To give an example, for the nucleon density 4p a = 0.06 fm~ 3 , a minimum of 
the energy expectation value (1731 was found at a = 0.616 fm _1 , b = 1.221, and 
c = —5.306, with a corresponding energy per a cluster of —9.763 MeV and a con- 
densate fraction of 0.750. The dependence of the condensate fraction on the nucleon 
density p = 4p a as determined in this exploratory calculation is displayed in Fig. [331 
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Fig. 33 Reduction of condensate fraction in a matter with increasing nucleon density. Exploratory 
calculations (full line) are compared with HNC calculations of Johnson and Clark 11231 (crosses). 
For comparison, we show estimates of the condensate fraction in the 0j (Hoyle) state of 12 C, 
according to Refs. |20 ;2TJ (stars). 
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Fig. 34 Occupation of the S orbital as a function of density using the3aOCMfor 12 C(2Tl. 



The reduction of the condensate fraction of a matter to roughly 0.8 as given by 
our calculation at nucleonic density 0.03 fm~ 3 agrees well with results of Suzuki 
EDI and Yamada ||2TI for 12 C in the Hoyle 0^ state. Using many-particle approaches 
to the ground-state wave function and to the THSR (Oj) state of 12 C, the occupa- 
tion of the inferred natural a orbitals is found to be quite different in the two cases. 
Roughly 1/3 shares (approaching equipartition) are found for the S, D, and G orbits 
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Fig. 35 Radial behaviors of the S orbit in the 12 C(0+) state with (a) R = 2.42 fm (p/p ~ 1.18), (b) 
R = 2.70 fm (p/po ~ 0.85), (c)R = 3.11 fm (p/p ~ 0.56), and (d) R = 4.84 fm (p/p - 0.15), 
where S denotes the nuclear radius of the 12 C(0 + ) state 1211 . 



in the ground (0j ) state, with a-cluster occupations of 1.07, 1.07, and 0.82, respec- 
tively (see Sec. |4.1.3"T l. On the other hand, in the Hoyle (0^) state, one sees enhanced 
occupation (2.38) of the S orbit and reduced occupation (0.29, 0. 16, respectively) of 
the D and G orbits (see also Sec. 14.1.3b . This corresponds to an enhancement of 
about 70% compared with equipartition. 

To get a more extended analysis, OCM calculations have been performed ||2TI 
for studying the density dependence of the S-orbit occupancy in the Hoyle state on 
the different densities p/po ~ CR(0j~ ) exp /R) 3 , in which the rms radius (R) of 12 C 
is taken as a parameter and /?(0^) exp =2.56 fm. A Pauli-principle respecting OCM 
basis f^ ) c J CM (v) with a size parameter v is used, in which the value of V is chosen 
to reproduce a given rms radius R of 12 C, and the a density matrix p(r,r') with 
respect to f^° CM (v) is diagonalized to obtain the 5-orbit occupancy in the + wave 
function. The results are shown in Fig. [34] The 5-orbit occupancy is 70 ~ 80 % 
around p/po ~ (^(0i h ) exp /^(0j) T HSR) 3 = 0-21, while it decreases with increasing 
p/po and amounts to about 30 ~ 40 % in the saturation density region. Figure [35] 
shows the radial behaviors of the S-orbit with given densities. A smooth transi- 
tion of the S-orbit is observed, with decreasing p / po, from a two-node 5-wave na- 
ture (p/po ~ 1.18) in Fig. l3"BT a) to the zero-node 5-wave one (p/po — 0.15) in 
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Fig. |35l d) ||2TI . The feature of the decrease of the enhanced occupation of the S 
orbit is in striking correspondence with the density dependence of the condensate 
fraction calculated for nuclear matter (see Fig. [33}. 

A more accurate and reliable variational description of a matter can be realized 
within the hypernetted-chain (HNC) approach to evaluate correlation integrals; this 
approach 112311261 largely overcomes the limitations of the cluster-expansion treat- 
ment, including the need for an explicit normalization constraint. Such an improved 
approach is certainly required near the saturation density of nuclear matter, where 
it predicts only a small condensate fraction 01231 . Of course, at high densities the 
simple Ali-Bodmer interaction ||93l ceases to be valid, and it becomes crucial to in- 
clude the effects of Pauli blocking. Once again, this conclusion reinforces the point 
of view that we can expect signatures of an a condensate only for dilute nuclei near 
the threshold of na decay. 



5.4 'Gap ' equation for quartet order parameter 

For macroscopic a condensation it is, of course, not conceivable to work with a 
number projected a particle condensate wave function as we did when in finite 
nuclei only a couple of a particles were present. We rather have to develop an anal- 
ogous procedure to BCS theory but generalized for quartets. In principle a wave 
function of the type \a) = exp[£i 2 34Zi234c|c2 f<: 3~ c 4 ]l vac ) would be the ideal gener- 
alization of the BCS wave function for the case of quartets. However, unfortunately, 
it is unknown so far (see, however, Ref. 01281 ) how to treat such a complicated 
many body wave function mathematically in a reasonable way. So, we rather at- 
tack the problem from the other end, that is with a Gorkov type of approach, well 
known from pairing but here extended to the quartet case. Since, naturally, the for- 
malism is complicated, we only will outline the main ideas and refer for details to 
the literature. 

Actually one part of the problem is written down easily. Let us guide from a 
particular form of the gap equation in the case of pairing. We have at zero tempera- 
ture II 10511 1061 

(£i+£ 2 )Kj 2 + (l -ni-n 2 )-Y,Vm'2'Ki>2' =2fiK n , (76) 

Z 1'2' 




Fig. 36 Graphic representation of the BCS mass operator in Eq. (1781 
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where K12 = {c\c%) is the pairing tensor, «,■ = (cfc t ) are the BCS occupation num- 
bers, and Vnvi 1 denotes the antisymmetrized matrix element of the two-body in- 
teraction. The £, are the usual mean field energies. Equation ((76) is equivalent to 
the usual gap equation in the case of zero total momentum and opposite spin, i.e. in 
short hand: 2=1 where the bar stands for 'time reversed conjugate'. The extension 
of d76l ) to the quartet case is formally written down without problem 

(£l234 -4^)K- 123 4 = (l-ni-M 2 )^£Vi21'2'K4'2'34 

1 1'2' 

+ (1 — ti\ — n-s)— ^I3i'3 /K "i'23'4 + all permutations. (77) 

2 j, 3 , 

with K4234 = (C1C2C3C4) the quartet order parameter. This is formally the same equa- 
tion as in Eq. d55l l with, however, the Fermi-Dirac occupation numbers replaced by 
the zero temperature quartet correlated single particle occupation numbers, similar 
to the BCS case. For the quartet case, the crux lies in the determination of those oc- 
cupation numbers. Let us again be guided by BCS theory or rather by the equivalent 
Gorkov approach 111291 . In the latter, there are two coupled equations, one for the 
normal single particle Green's function (GF) and the other for the anomalous GF. 
Eliminating the one for the anomalous GF in inserting it into the first equation leads 
to a Dyson equation with a single particle mass operator, 

M?.5 S (©)=£^^ with A n = -Wvn,34M. (78) 

This can be graphically represented in Fig. [36] where (cc) stands for the order pa- 
rameter K\2 and the dot for the two body interaction. 

The generalization to the quartet case is considerably more complicated but 
schematically the corresponding mass operator in the single particle Dyson equa- 
tion can be represented graphically as in Fig. [37] with the quartet order parameter 
(cccc). Put aside the difficulty to derive a manageable expression for this 'quartet' 
single-particle mass operator, what immediately strikes is that instead of only one 
'backward going line' with (— p,— a) as in the pairing case, we now have three 
backwards going lines. As a consequence, the three momenta ki, k2, k3 in these 
lines are only constrained so that their sum be equal to kj + k2 + k3 = — p and, 
thus, the remaining freedom has to be summed over. This is in strong contrast to the 
pairing case where the single backward going line is constrained by momentum con- 
servation to —p. So, no internal summation occurs in the mass operator belonging 
to pairing. The consequence of this additional momentum summation in the mass 
operator for quartetting leads with respect to pairing to a completely different ana- 
lytic structure of the mass operator in case of quartetting. This is best studied with 
the so-called three hole level density #3/, (ft)) which is related to the imaginary part 
of the three hole Green' s function G 3h (k u k 2l k i ;(o) = (/i/ 2 / 3 + /1/2/3 )/(<» + £123 ) 
with £ 12 3 = £1 + £2 + £3 and / = 1 — / by 



Nuclear Alpha-Particle Condensates 



61 




. , f d 3 k\ d 3 k2 d 3 k$ , , 

83h{co) = -J (my (my (2n/ mG ] ^ k ^+^) 



I 



d^h d 3 k 2 d 3 k 3 



(2tt) 3 (2tt) 3 {2nf 
x {fihh+fihh)nS{<d + ei + e 2 + £3).) (79) 

In Fig.[38]we show the level density at zero temperature (f(co) = 9(—0))), where 
it is calculated with the proton mass m = 938.27 MeV (natural units) (52j. Two cases 
have to be considered, chemical potential jj. positive or negative. In the latter case we 
have binding of the quartet. Let us first discuss the case ji>0. We remark that in this 
case, the 3h level density goes through zero at CO = 0, i.e., since we are measuring 
energies with respect to the chemical potential ji, just in the region where the quartet 
correlations should appear. This is a strong difference with the pairing case where 
the Hi level density, g lh (co) = J (f k + f k )8((Q + S k ) = J jj^8{co + e k ), does 
not feel any influence from the medium and, therefore, the corresponding level den- 
sity varies (neglecting the mean field for the sake of the argument) like in free space 
with the square root of energy. In particular, this means that the level density is finite 
at the Fermi level. This is a dramatic difference with the quartet case and explains 
why Cooper pairs can strongly overlap whereas for quartets this is impossible as we 
will see below. We also would like to point out that the 3/j level density is just the 
mirror to the 3p level density which has been discussed in Ref. 11301 . 

For the case where ji < there is nothing very special, besides the fact that it 
only is non- vanishing for negative values of (0 and that the upper boundary is given 
by co = 3jJ.. Therefore, the level density of Eq. ((79) is zero for co > 3fl. 

With these preliminary but crucial considerations we now pass to the evaluation 
of the single-particle mass operator with quartet condensation. Its expression can be 
shown to be of the following form 

quartet, s = y ^mAjflhh+flhfAj^A ( g Q) 
U 234 » + £234 

with 

^1234 = 2^12,l'2'533'5 4 4/(c 1 /C2'C 3 /C4/). (81) 
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Fig. 38 3h level densities defined in Eq. i79\ for various values of the chemical potential fl at zero 
temperature 1521 . 



Again, comparing the quartet single-particle mass operator (f80b with the pairing 
one d78l ), we notice the presence of the phase space factors in the former case while 
in Eq. d78l l they are absent. As already indicated above, this fact implies in the 
quartet case that only the Bose-Einstein condensation phase is born out whereas a 
'BCS phase' (long coherence length) is absent. The complexity of the calculation in 
Eq. (f80b is much reduced using for the order parameter (cccc) our mean field ansatz 
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Fig. 39 Single particle wave function <p(k) in £-space (left), for r-space <p(r) (middle), and 
occupation numbers (right) at jit = —5.26 (top), —1.63 (middle), and 0.55 (bottom), with zero 
temperature. The r-space wave function (f>(r) is derived from the Fourier transform of <p(k) by 
(f>(r) = J d 3 ke' k ' r (p(k)/(27z) 3 . The dashed line in the left panels correspond to the Gaussian with 
same norm and r.m.s. momentum as (p(k) 1521 



projected on zero total momentum, as it was already very successfully employed 
with Eq. ([58), 

^k 1 k 2 ,k 3 k 4 = <KN)<KN)<p(|k 3 |)<p(|k 4 |) 

x(27r) 3 5(k 1 +k2 + k 3 + k 4 ), (82) 
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Fig. 40 -ImMl uar,et (/ti , ffl + irj) in Eq. {80j as a function of CO for fl = -4.9MeV (left) and for 
fl = 0.55MeV (right) at zero temperature |52| . 



where Xo is the spin-isospin singlet wave function. It should be pointed out that 
this product ansatz with four identical OS single particle wave functions is typical 
for a ground state configuration of the a particle. Excited configurations with wave 
functions of higher nodal structures may eventually be envisaged for other physical 
situations. We also would like to mention that the momentum conserving 8 function 
induces strong correlations among the four particles and (l82l is, therefore, a rather 
non trivial variational wave function. 

For the two-body interaction of V\2 V2' m Eq. d8TT l, we employ the same separable 
form (59) as done already for the quartet critical temperature. 

At first let us mention that in this pilot application of our selfconsistent quartet 
theory, we only will consider the zero temperature case. As a definite physical ex- 
ample, we will treat the case of nuclear physics with the particularly strongly bound 
quartet, the a particle. It should be pointed out, however, that if scaled appropriately 
all energies and lengths can be transformed to other physical systems. For the nu- 
clear case it is convenient to measure energies in Fermi energies £f — 35 MeV and 
lengths in inverse Fermi momentum kj l = 1.35 -1 fm. 

The single particle wave functions and occupation numbers obtained from the 
above cycle are shown in Fig. [39] We also insert the Gaussian wave function with 
same r.m.s. momentum as the single particle wave function in the left figures in 
Fig. [39] As shown in Fig. [39] the single particle wave function is sharper than a 
Gaussian. 

We could not obtain a convergent solution for fl > 0.55 MeV. This difficulty is 
of the same origin as in the case of our calculation of the critical temperature for a 
particle condensation. In the r.h.s. panels of Fig. [39] we also show the corresponding 
occupation numbers. We see that they are very small. However, they increase for 
increasing values of the chemical potential. For fl = 0.55 MeV the maximum of the 
occupation still only attains 0.35 what is far away from the saturation value of one. 
What really happens for larger values of the chemical potential, is unclear. Surely, 
as discussed in Sec. 15.21 the situation for the quartet case is completely different 
from the standard pairing case. This is due to the fact, as already mentioned, that the 
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Fig. 41 -2lmG(k, co + irj) in Eq. l80t as function of CO for fl 
0.55 MeV (bottom) at a zero temperature 1521 . 
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3h level density goes through zero at 0) = 0, i.e. just at the place where the quartet 
correlation should build up for positive values of jj.. Due to this fact, the inhibition 
to go into the positive /i regime is here even stronger than in the case of the critical 
temperature ||50) . 

The situation in the quartet case is also in so far much different, as the 3h Green's 
function produces a considerable imaginary part of the mass operator. Figure |4Q] 
shows the imaginary part of the approximate quartet mass operator of Eq. d80b for 
H < and n > 0. These large values of the damping rate imply a strong violation of 
the quasiparticle picture. In Fig.|4T]we show the spectral function of the single par- 
ticle GF. Contrary to the pairing case with its sharp quasiparticle pole, we here only 
find a very broad distribution, implying that the quasiparticle picture is completely 
destroyed. How to formulate a theory which goes continuously from the quartet case 
into the pairing case, is an open question. One solution could be to start right from 
the beginning with an in medium four body equation which contains a superfiuid 
phase. When the quartet phase disappears, the superfiuid phase may remain. Such 
investigations shall be done in the future. 
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6 Summary and conclusions 

We discussed a condensation in nuclear systems. One remarkable manifestation is 
the Hoyle state (Oj) in 12 C at 7.65 MeV with a gas-like structure of three a-particles, 
trapped by a shallow self-consistent mean field of wide extension, in which the cm. 
motion of the a particles occurs dominantly in the lowest OS orbit. We found that 
a simple wave function of the a-condensate type, called the THSR wave function, 
describes very nicely the structure of the Hoyle state and reproduces the inelastic 
form factor of 12 C(e,e') and others quantities. The condensate feature of the Hoyle 
state was confirmed by the calculation of the bosonic occupation numbers in diag- 
onalizing the bosonic density matrix. It was shown that the occupation of the OS 
state of the a-particles is over 70 % for the Hoyle state, and the remaining compo- 
nent (30 %) comes from residual correlations, mostly of the Pauli type, among the 
a particles. In spite of the very different number of particles and other important 
differences, the situation has some analogy with the case of cold bosonic atoms. 

We conjectured that the a-particle condensates also exist in heavier self-conjugate 
nuclei. Theoretical calculations of the OCM type indicate that the 0g state at 15.1 
MeV in ie O is a strong candidate. So far we do not dispose of sufficient experi- 
mental data to confirm its nature. Experiments are under way and being analyzed. 
This analogue of the Hoyle state in ie O has many similarities with the original one 
in 12 C: it lies a couple of hundred keV above the 4a disintegration threshold. It is 
quite strongly excited by (e,e'). Its width is 160 keV. This is much larger than for the 
Hoyle state in 12 C but with respect to its energy it is still unusually small. The large 
width stems from a position higher up in Coulomb barrier and also the Coulomb 
barrier itself has become slightly lower. The situation in ie O with respect to alpha 
clustering is considerably more complicated than in 12 C. Results from the 4a OCM 
calculations showed that 2nd up to 5-th + states in ie O have a+ 12 C structures. 
Only the 6-th + state is the analogue to Hoyle state. We also discussed the results 
of the THSR wave function for le O. 

As for the heavier a-particle condensates, we found first that they are predicted 
to be slightly above their na threshold in the A = 4« nuclei but below the Coulomb 
barrier, and second the phenomenon will terminate at about eight to ten a's as the 
confining Coulomb barrier fades away. However, the concept of a condensation in 
nuclei can be generalized to non self-conjugated nuclei (A ^ 4n). Since the nuclear 
a-particle condensation is described dominantly as a product state of a particles 
occupying the lowest OS orbit, the counterpart in A ^ An nuclei should still be pre- 
sented as a product state of the constituent clusters in the 05 state. For instance, 
we can conjecture product states composed of a's, a few neutrons and/or s-wave 
clusters (d, t, 3 He) such as (OS) 2 (OS), in n B and (0S) 3 a (0S) n in 13 C etc. Indeed, 
our OCM calculations indicate that they appear slightly above their three- and four- 
cluster disintegrated thresholds, 2a + t for n B and 3a + n for 13 C, as positive -parity 
states with J % = l/2 + . These results encourage us to conjecture that cluster-gas -like 
states described by antisymmetrized product wave functions of constituent clusters, 
all in the OS level, can exist in general in excited states of low density in light nuclei. 
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We dwell on the fact that concepts developed for infinite nuclear matter are of 
value also to interpret properties in finite nuclei and to construct useful approxi- 
mations. As examples, we refer to pairing, two and more body correlations, and 
one body occupation numbers. Pairing definitely also is a useful concept for many 
finite nuclei, in spite of the fact that nuclei are by far not macroscopic objects. 
For example, the strong reduction of measured moments of inertia of such nuclei 
compared with the classical values are explained as a consequence of superfluid- 
ity |[T9lll31l[l32l . In this sense, we discussed nuclear a-particle condensation as 
the analogue to the number-projected BCS wave function, replacing Cooper pairs 
by a particles. A real macroscopic phase of condensed a's may be formed during 
the cooling process of compact stars 111331 . where one predicts the presence of a- 
particle condensates 01211 . On the other hand, a possibility of quartetting with cold 
atoms in which fermions are trapped in four different magnetic substates also have 
been discussed 01 13111 15 134|. Theoretical and experimental works in this direction 
will also be useful and helpful to investigate the low-density bosonic a-particle gas 
states in nuclei. 

In conclusion, the idea of a-particle condensation in nuclei is novel. A com- 
pletely new nuclear phase in which a particles move like in a gas as quasi- 
elementary constituents is surely intriguing. In order to bring deeper insights into 
the role of clustering and quantum condensates in the systems of strongly inter- 
acting fermions, it is hoped that more a-particle states in nuclei and/or many a's 
around a nuclear core, including cluster-gas-like states composed of a's, f's and n's 
etc., will be observed in the near future. 
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